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SUMMARY 


The Extended Performance Technology study was an 
investigation of advanced discharge chambers and thruster 
components that were designed to operate under conditions that 
result in an increase in the thrust and thrust-to-power ratio of 
the state-of-the-art J-series thruster. Achieving the extended- 
performance capability required a thruster design that would 
operate efficiently at beam currents as high as = 5 A and at 
net-to-total accelerating-voltage ratios as low as R = 0.2. The 
performance goals were achieved using a ring-cusp-magnetic- 
confinement thruster design that demonstrated a baseline ion- 
production cost of only = 125 eV/ion, a beam current as high 
as J 5 = 5.3 A, and an acceptable range of net-to-total 
accelerating- voltage ratios of 0.2 < R < 0.8. When compared with 
the performance of the state-of-the-art J-series thruster, these 
performance levels reoresent nearly a f actor-of-three increase in 
thrust, a factor-of-two increase in thrust-to-power ratio, and a 
35% reduction in baseline ion-oroduction cost. 

The extended operating range and performance objectives 
translate into exacting requirements on the discharge chamber, 
ion-extraction assembly, hollow cathode, and propellant-flow 
electrical isolators. Therefore, the program was logically 
divided into four major elements or tasks that dealt with these 
critical technology areas. 

In the first element, which consisted of three related 
tasks, several advanced discharge— chamber configurations and 
their variations were investigated. Each of the designs 
demonstrated performance that was comparable to that of the 
highly optimized J-series thruster, but the ring-cusp design was 
found to be far superior over the other two (a single-cusp and a 
multipole magnetic-confinement geometry). Extensive performance 
and plasma-diagnostic measurements were conducted using a 
laboratory model of the ring-cusp thruster, leading to a 
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preliminary description or model of the discharge plasma that is 
consistent with its high-performance characteristics. Our 
plasma-property measurements indicate that the ring-cusp- 
magnet ic-conf inement geometry effectively concentrates the 
ionization region near the ion-extraction assembly, resulting in 
a high probability that the ions produced in this volume 
will be extracted into the ion beam. 

The investigation of a flexible-mount electrode-supporting 
structure that was developed under a previous NASA contract (NAS 
3-21040), demonstrated that an advanced three-grid ion-extraction 
assembly could operate under high-beam-current extraction 
conditions (high power), while at the same time maintaining the 
ratio of net-to-total accelerating voltage as low as R = 0.2. 

The technology of three-grid- ion-extraction assemblies was 
further developed under this contract by performing a parametric 
analysis of the electrode-aperture displacement required to 
vector individual beamlets in a direction to minimize off-axis 
thrust losses. 

A model of hollow-cathode lifetime that was also developed 
under NASA contract NAS 3-21040 was used to design a cathode that 
can operate under high-emission-current conditions for extended 
periods of time. Over 60,000 A-hr of operation were demonstrated 
with scaled versions of the high-performance cathode used in the 
state-of-the-art J-series thruster; our test results indicate 
that long lifetime can be expected at emission currents as high 
as J E * 40 A. This will be possible as long as the diameter of 
the flow-restricting orifice is large enough to maintain the 
ratio of emission current to aperture diameter at about the same 
value as that of the J-series design. 

A study of the breakdown- volt age characteristics of the 
state-of-the-art propellant- flow electrical isolator established 
the maximum conservative isolation voltage as 340 V per section. 
Our measurements further showed that the total breakdown voltage 
is directly proportional to the number of individual isolator 
sections. This finding enables a specified voltage standoff to 
be achieved simply by scaling the existing isolator design so 
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that the product of the number of sections and the section rating 
of 340 V provides a comfortable margin over the required value. 

Other tasks included the design, fabrication, and testing of 
a simplified vaporizer power supply and a comprehensive 
investigation of the uncertainties in measuring thrust-loss 
factors. The power-supply simplifications incorporated 
multiplexing techniques to eliminate over 30% of the components 
required in the three vaporizer power supplies used in the power 
orocessor for the J-series thruster. The results of the 
uncertainty analysis indicated that the overall error in the 
thrust-loss factors, a and F^, is less than 1.4%. 


FOREWORD 


The work described in this report was performed during a 
period of approximately three years at Hughes Research 
Laboratories (HRL). The technical effort was initially 
conducted in the High Voltage Technology Department, managed by 
Dr. H.J. King. In May 1981 this department was reorganized and 
became the Plasma Physics Department, managed by Dr. J. Hyman. 
The program was managed by Dr. J.R. Beattie of HRL and 
monitored initially by Mr. W.R. Kerslake, and later by 
Mr. V.K. Rawlin, both of NASA-Lewis Research Center (LeRC). 
Major technical contributions to this effort were made by: 
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SECTION 1 


INTRODUCTION 

By the beginning of the present program the 30-cm-diameter 
J-series thruster had reached a level of development that would 
enable confident application of prime-propulsion technology, and 
advanced studies had shown promise for extending the performance 
of the baseline design. The objective of this program was, 
therefore, to investigate discharge-chamber, ion-accelerator, 
and component designs that would produce significantly higher 
thrust and higher thrust-to-power ratio than could be achieved 
with the baseline 30-cm thruster. 

A. PROGRAM GOALS 

The work was performed under seven technical tasks. Tasks 
I, VIII, and IX required evaluation of advanced discharge- 
chamber and ion-extraction assemblies that would enable 
operation at increased thrust and increased thrust— to— power 
ratio. Task II required investigation of simplifications that 
would reduce the complexity and cost of the power processors 
required for the extended-performance thruster. Task III 
required lifetime evaluation of a high-emission-current 
(j E = 40 a) hollow cathode. Task IV required measurements of 
the pressure-dependent electrical-breakdown voltage of 
individual sections of the state-of-the-art propellant-flow 
electrical isolator. Task X required calculations of the 
electrode compensation necessary to minimize the off-axis thrust 
loss of three-grid ion-extraction assemblies. Other work 
required under Task I included a detailed description of the 
technique used for determining the thrust-loss factors o and 
F t , and an estimate of the uncertainties in the measurements. 
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B. PROGRAM ACCOMPLISHMENTS 


Under Tasks I , VIII, and IX, three discharge chambers and 
two ion-extraction-assembly designs were investigated. The 
discharge-chamber designs included the ring-cusp, single-cusp, 
and multipole magnetic-confinement arrangements. The ion- 
extraction assemblies included a three-electrode arrangement 
that employs flexible mounts for supporting each grid and a two- 
electrode design that has the same aperture dimensions as the J- 
series ion-extraction assembly, but which uses flexible mounts 
for supporting the accelerator electrode. The performance of 
these discharge-chamber and ion-extraction designs was 
documented under the program. The ring-cusp configuration was 
found to give the best discharge-chamber performance, having a 
baseline ion-production cost of about e* = 125 eV/ion. The 
two-grid and three-grid ion-extraction assemblies were 
successfully operated at beam currents as high as = 5 a, 
and the three-grid ion-extraction assembly was operated at net- 
to-total accelerating-voltage ratios in the range of 0.1 < r < 
0.8. A comprehensive analytical and experimental investigation 
into the uncertainties in the measured values of the thurst-loss 
factors o and F^. was completed, and the results showed that 
the uncertainty in these factors is less than 1.5%. 

Under Task II, several power-processor simplifications were 
investigated. Of the various design simplifications that we 
identified, a multiplexed vaporizer power supply was selected 
for fabrication and testing with a J-series-equivalent thruster. 
The time-sharing device performed very well, and the simplified 
design results in a significant reduction (= 70%) in parts count 
(and, therefore, lower cost and higher reliability) when 
compared with the three individual vaporizer power supplies 
required for the J-seriea thruster. 
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Work under Task III included the design of a high-emission- 
current (J E = 40 A) hollow cathode. Three of these cathodes 
were tested for 500 hr each, accumulating 20,000 A-hr of 
operation per cathode. The lifetest results showed that the 
aperture wear is excessive unless the diameter of the flow- 
restricting orifice is d Q > 1.27 mm. Temperature measurements 
showed that the lifetime of the high-emission-current cathode 
(based on depletion of barium from the porous-tungsten insert) 
should be comparable to that of the J-series insert ( “30,000 hr). 

Under Task IV, the breakdown voltages of individual 
sections of a propellant-flow electrical isolator were measured 
as a function of flow rate and isolator-base temperature. The 
minimum breakdown voltage was found to be about 340 V per 
section, and the total breakdown voltage was found to scale 
linearly with the number of isolator sections. 

Work under Task X resulted in design curves for use in 
determining the electrode-aperture displacement required to 
vector the beamlets emerging from three-grid ion-extraction 
assemblies. The parametric results cover a wide range of design 
variables, including electrode thickness, aperture diameter, and 
interelectrode spacing. 

The work performed under the contract is described in 
detail in the following sections. The report also describes 
plasma-diagnostic measurements and analyses that were performed 
under Hughes funding, enabling the operation and plasma physics 
of the ring-cusp discharge chamber to be more fully documented. 
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DISCHARGE-CHAMBER INVESTIGATION 
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THRUST-TO-POWER RATIO, F/P in mN/kW 


o 


SPECIFIC IMPULSE, l sp , sec 


Figure 2-1 


Variation of thrust- to- power ratio with 
specific impulse for an electrostatic 
ion thruster. Calculations represent 
the idealized case, where fixed-power, 
thrust, and un-ionized propellant losses 
have been neglected. 





INCREASE INTHRUST-TO-POWER RATIO, A(F/P in ), % 


ej = 100 eV/ion 



NET-TO-TOTAL ACCELERATING - VOLTAGE RATIO, R 


Figure 2-2. 


Performance gains that could be achieved by 
reducing the specific discharge power, e^* 
and net-to-total accelerating-voltage ratio 
R, of the state-of-the-art 30-cm thruster. 

In performing these calculations, we have 
neglected fixed-power losses (vaporizers, 
keepers, etc.) and we have used the thrust- 


loss factor* y, propellant-utilization 
efficiency, n m , and power processor 


efficiency, npp indicated in the figure. 




gain is expressed as an increase in thrust-to-power ratio over 
that of the J-series thruster. Figure 2-2 indicates that a 60% 
increase in the performance of the baseline thruster can be 
achieved by lowering the ion-production cost to = 135 eV/ion 
and by reducing the net-to-total accelerating-voltage ratio to 
R = 0.2. 

Achieving the desired performance improvements of increased 
thrust and increased thrust-to-power ratio required a discharge 
chamber that could produce a sufficiently dense plasma with mini- 
mum power consumption, and an ion-extraction assembly that could 
efficiently extract the beam current at minimum net-to-total 
accelerating-voltage ratio. These requirements were demonstrated 
using the advanced-technology ring-cusp discharge chamber des- 
cribed below, and the advanced-technology three-grid ion- 
exJtraction assembly described in Section 3. 

A. RING-CUSP DISCHARGE CHAMBER 

Figure 2-3 illustrates the high-performance, laboratory-model 
ring-cusp discharge chamber that we fabricated and tested under 
this program. The design is based on the ring-cusp magnetic- 
confinement configuration developed at NASA-Lewis Research 
Center. 2-3 The entire cylindrical sidewall and circular endwall 
are maintained at anode potential and were fabricated from 
commercially available iron pipe and sheet stock. The cathode and 
cathode magnet are electrically isolated from each other and from 
the endwall and (in the laboratory-model thruster) can be 
independently moved during thruster operation using the cathode 
and magnet gear-drive arrangement that we also designed and 
fabricated.* The main propellant flow is introduced into the 
discharge chamber through a plenum which is electrically isolated 
from the endwall. A photograph showing the interior of the 
laboratory-model ring-cusp discharge chamber is presented in 
Figure 2-4. 


*Details of the gear-drive arrangement are provided in 
Section 2-B. 
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Piqure 2-4. Photograph of the laboratory-model 
ring-cusp discharge chamber. 



High-f leld-strength samarium cobalt (SmCo 5 ) magnets are used 
to produce the cusp fields along the boundary, and an array of 
Alnico-V magnets surrounds the cathode. The resultant confining 
magnetic field has the divergent-cusp shape evident in the iron- 
filings map of Figure 2-5. The SmCo 5 magnets are held in place by 
the strong forces of attraction between them and the soft-iron 
boundary, the magnetic-field strength at the surface of the 
magnets is approximately 2,400 G. 


The electrically isolated magnet covers (shown in the photo- 
graph of Figure 2-4 and the schematic of Figure 2-6) allow the 
anode-current distribution to be determined, and they provide 
thermal protection for the magnets which have relatively low 
electrical conductivity and a maximum temperature rating of 325®c, 
A black-oxide coating on the exterior of the discharge chamber 
also aids in keeping the magnets cool by enhancing the radiative 
transfer to the surroundings, with this arrangement, steady-state 
testing at discharge powers up to 750 w is possible, while main- 
taining the magnets -65*c below their maximum temperature rating 

We measured the radial and axial components of the magnetic 
field produced by the "rings" of SmCo s magnets and the cathode 
magnet, using the 2-axis gaussmeter probe and positioning mecha- 
nism shown in Figure 2-7. The resultant magnetic field was calcu- 
ated from the measured components of B, and contours of the 
scalar fields are shown in Figures 2-5 and 2-8 to lie approxi- 
mately normal to the field lines obtained from iron filings and 
parallel to the chamber boundaries. 

The axial magnetic field on the discharge-chamber axis can be 
varied during thruster operation using the movable cathode-magnet 
arrangement shown in Figure 2-3. The axial field is produced by 
the contributions of the SmCo 5 magnets located on the chamber 
sidewall and endwall, and the cathode magnet, which is a cylindri- 
cal array of Alnico-V rods with an iron end-cap used to concen- 
-rate the field on axis. The variation of the axial magnetic 
field on the chamber axis is presented in Figure 2-9, with the 
location (L) of the movable cathode magnet as a parameter. 
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Figure 2-8. Computer-generated magentic-f ield contours 
obtained from magnetic-field measurements 
performed in the ring-cusp discharge chamber. 





> J • ariv 




33 







Figure 2-10 presents a schematic that shows the power- 
supply arrangement used in testing the thruster and defines the 
electrical quantities that were measured. Equations used in 
calculating thruster performance are summarized in Table 2-1. 

1 . Performance Characteristics 

Initial testing of the ring-cusp discharge chamber 
established the minimum-emission-current requirement (minimum 
eV/ion) for discharge voltages in the range 28 V < Vq < 44 V. 
Figure 2-11 presents the variation of cathode- and main-vaporizer 
temperatures* with emission current, Je, and discharge 
voltage, Vq, for a beam current of Jb ■ 2 A. These 
measurements documented the current-voltage characteristics of the 
discharge chamber, and they established the sensitivity of the 
cathode and main flow rates to the operating parameters, Je and 
Vq. The "knees" of the curves define the minimum emission 
current (minimum e^), and from Figure 2—11 it is evident that 
they also correspond to a rapid increase in main propellant flow 
and a rapid decrease in cathode flow. The baseline-performance 
level is also apparent in the variation of accelerator current, 
j A , with emission current, as shown in Figure 2-12. The 
vertical part of these curves defines the minimum ion-production 
cost, which lies in the range 120 < ei < 154 eV/ion for 
discharge voltages in the range 44 V > Wq > 28 V. The 
performance improvement at higher discharge voltages is attributed 
to the larger ionization cross sections at the higher electron 
energies and more-efficient ion production associated with the 
reduction in cathode flow rate (that accompanies operation at 
higher discharge voltage). 


*After obtaining these results, it was determined that 
propellant was leaking from the base of the main vaporizer. 
Therefore, the results of Figures 2-11 and Figures 2-13 to 2-15 
give correct trends, but they are not absolute measurements. 
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Table 2-1. Summary of Equations used to Calculate 
Performance of Ring-Cusp Thruster 
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Figure 2-11. Variation of cathode and main— vaporizer 
temperature with cathode-emission cur- 
rent (Jk = 2 A) . 
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2. Optimization of the Cathode/Keeper Assembly and 

Cathode-Magnet Locations 

The effects of the locations of the cathode/keeper assembly 
and cathode magnet are illustrated in Figure 2-13, which shows 
cathode- and main-vaporizer temperatures as a function of the 
axial position of these components. The temperature measurements 
show that for a given location of the cathode-keeper assembly, 
the magnet position has little, if any, effect on main flow rate, 
while the cathode flow rate is minimized by increasing the dis- 
tance between the cathode magnet and cathode orifice. The B 
measurements presented in Figure 2-9 and the temperature/position 
curves of Figure 2-13 allow one to infer that lower cathode flows 
are associated with a magnetic-field gradient that points toward 
the cathode, while the higher flows are associated with a gradient 
that points away from the cathode. Because of the positive and 
negative slopes of the vaporizer-temperature curves of Figure 2- 
13, it is natural to question whether or not an optimum location 
of the cathode/keeper and magnet exists which would minimize the 
sum of the cathode- and main-vaporizer flow rates, maximizing the 
propellant-utilization efficiency. Using the vaporizer- 
calibration curves of Figure 2-14, we converted the temperature 
data of Figure 2-13 to flow rates and then used these calibrated 
flows to calculate propellant-utilization efficiencies. Our 
results are presented in Figure 2-15, which shows that the 
relative propellant-utilization efficiency is a maximum for the 
cathode/keeper located nearest the endwall. With the location of 
the cathode defining the "virtual" endwall, or chamber length, 
this finding is equivalent to achieving increased propellent 
utilization with increased neutral residence time. 

3. Extended-Performance Operation 

Having demonstrated superior performance (low baseline 
e^) of the ring-cusp thruster design over that of the state- 
of-the-art J-series thruster for operation at the nominal design 


APORlZER TEMPERATURE, °C 




CATHODE 

VAPORIZER 


V D = 33 V 
J E *9.5 A 
Jb = 2 A 


MAGNET 

~ SYMBOL POSITION, cm 


♦X ♦ 


MAIN VAPORIZER 


CATHODE POSITION, cm 


Figure 2-13. Summary of cathode and main-vaporizer- 
temperature measurements obtained using 
the movable cathode/cathode-magnet 
arrangement in the ring-cusp thruster. 
(Dimensions are relative to the 
thruster endwall). 






CATHODE FLOW RATE, M CATH , A MAIN FLOW RATE, M MA j N , 


1.70 1.72 1.74 1.76 1.78 1.80 1.82 1.84 x 10~ 3 

RECIPROCAL OF MAIN-VAPORIZER TEMPERATURE, T mV -1 , °K~ 1 
a. MAIN VAPORIZER 


1.58 1.60 1.62 1.64 1.66 1.68 1.70 1.72x10 -3 

RECIPROCAL OF CATHODE-VAPORIZER TEMPERATURE, T^" 1 , 0 K -1 
b. CATHODE VAPORIZER 


Figure 2-14. Vaporizer calibration curves. 

(Not absolute measurements ) 
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point of Jjj - 2 A, we investigated its performance at higher 
power levels. Figure 2-16 presents performance measurements for 
the ring-cusp thruster operated over the beam-current range 
2 A < J b < 5 A. These results suggest a "knee" performance of 
about 125 eV/ ion and show that little gain in propellant 
utilization occurs above 180-200 eV/ion. The increase in 
maximum-propellant-utilization efficiency with beam current is 
consistent with Kaufman's theory of constant neutral loss rate. 
The 4-A and 5-A curves do not cover the range > 180 eV/ion 
because the maximum power that the laboratory-model thruster can 
dissipate without exceeding 260°C magnet temperature is about 
750 w. This arbitrary upper limit on allowable magnet 
temperature was selected to provide a comfortable margin-of- 
safety below the maximum safe value of 325°C. 

The propellant-utilization efficiencies of Figure 2-16 have 
not been corrected for the effects of multiply charged ions, so 
that values of n m greater than 100% are indicated. Inasmuch 
as the effects of multiply charged ions are expected to be 
relatively small at lower discharge voltages and beam currents, 
the high propellant-utilization efficiencies observed at Vjj = 

32 V and J b = 2 A suggest that the ring-cusp discharge chamber 
is also highly efficient in this aspect of its performance. 

Table 2-2 Summarizes the electrical measurements and 
performance calculated from these measurements for operation 
near the "knee" of the performance curves. The (uncorrected) 
thrust corresponding to a beam current of J b ■ 5 A is about a 
factor-of-three greater than that produced by the J-series 
thruster. 

Beam-Current-Dens itv Profile 

We operated the ring-cusp discharge chamber with a three- 
grid ion-extraction assembly (described in the Section 3) for 
beam currents up to J b = 5.3 A. The extracted ion-current- 
density profiles measured 3-mm downstream of the decelerator 
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Figure 2-16. Ring-cusp-thruster performance measurements 
(Propellant-utilization efficiencies 
uncorrected for multiply charged ions.) 
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Table 2-2. Summary of Electrical and Performance 

Measurements for the Ring-Cusp Thruster 



TEST NUMBER 



v b 

V 

1200 

1200 

1200 

1200 

1325 

1325 

1325 

1290 

1290 

1290 

1225 

1225 

J b 

A 

2.0 

2.0 

2.0 

2.0 

3.0 

3.0 

3.0 

4.0 

4.0 

4.0 

5.0 

5.0 

V D 

V 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

36.0 

J D 

A 

9.0 

10.0 

12.0 

14.0 

15.0 

18.0 

21.0 

20.0 

22.0 

24.0 

25.0 

27.0 

J E 

A 

7.0 

8.0 

10.0 

12.0 

12.0 

15.0 

18.0 

16.0 

18.0 

20.0 

20.0 

22.0 

J DRCEL 

mA 

4.95 

1.17 

1.19 

1.50 

1.98 

2.13 

2.33 

3.13 

3.18 

3.50 

5.40 

4.47 

V CK 

V 

13.0 

10.9 

10.0 

9.44 

10.2 

9.69 

8.90 

9.96 

9.29 

9.20 

9.70 

9.23 

J CK 

A 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

V 

V ACCEL 

V 

-300 

-300 

-300 

-300 

-525 

-525 

-523 

-860 

-860 

-860 

-1225 

-1125 

J ACC.3L 

mA 

14.0 

4.83 

4.08 

3.88 

8.35 

6.72 

5.74 

13.75 

10.6 

8.98 

21.5 

13.1 

V NK 

V 

15.5 

15.5 

15.5 

15.5 

15.3 

15.3 

15.5 

14.0 

13.9 

13.9 

14.4 

14.5 

J NK 

A 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.5 

2.5 

2.5 

3.1 

3.1 

V G 

V 

11.0 

10.3 

10.2 

10.2 

12.4 

12.2 

12.4 

11.3 

12.3 

12.3 

12.7 

13.3 

"W 

eq A 

2.894 

1.818 

1.629 

1.524 

2.836 

2.452 

2.286 

4.122 

3.392 

3.192 

4.361 

4.167 

"tv 

eq A 

0.090 

0.139 

0.212 

0.30C, 

0.164 

0.271 

0.403 

0.182 

0.265 

0.341 

0.207 

0.286 

"tOTAL 

eq A 

2.984 

1.957 

1.841 

1.832 

3.000 

2.723 

2.685 

4.304 

3.657 

3.533 

4.568 

4.453 

n mD (UNC) 

% 

67.0 

102.2 

108.6 

109.2 

100.0 

110.2 

111.6 

92.9 

109.4 

113.2 

109.5 

112.3 

P b 

W 

2400 

2400 

2400 

2400 

3975 

3975 

3975 

5160 

5160 

5160 

6125 

6125 

P (1) 
TOT 

W 

2718 

2751 

2821 

2893 

4485 

4592 

4700 

5826 

5901 

5973 

6963 

7038 

n e (D 

% 

88.3 

87.3 

85.1 

83.0 

88.6 

86.6 

84.6 

88.6 

87.4 

86.4 

88.0 

87.0 

e i 

W/A 

126 

144 

180 

216 

144 

180 

216 

144 

162 

180 

144 

158 

p(2) 

mN 

141.3 

141.3 

141.3 

141.3 

222.7 

222.7 

222.7 

292.9 

292.9 

292.9 

356.8 

356.8 

T (2,3) 
SP 

sec 

2323 

3544 

3765 

3786 

3643 

4015 

4066 

3340 

3933 

4069 

3836 

3934 

n (1,2,3) 
n T 

% 

59.2 

89.2 

92.4 

90.6 

88.6 

95.4 

94.4 

02.3 

95.6 

97.8 

96.4 

97.7 


1. Neglects vaporizer power 

2. oF t « 1.0 

3. Neglects neutralizer flow 


A, 














electrode are presented in Figure 2-17 for beam currents in the 
ranqe 1 A < < 5 A. The calculated beam-flatness parameter 

for these profiles is approximately F = 0.4, which is about 
equal to that of divergent-field designs such as the 30-cm J- 
series thruster. The relatively non-uniform profile (as com- 
pared to those of the single-cusp and multipole discharge cham- 
bers discussed below) is believed to be a result of the strong 
magnetic field, which penetrates into the bulk of the discharge- 
chamber volume, giving rise to plasma-density gradients. 

5 . Wall-Frobe Measurements 

We measured the distribution of electron and ion currents 
to the chamber boundaries using the anode-cover arrangement of 
Figure 2-6 to measure the electron currents to the magnet rings, 
and using the ion-flux probe illustrated in Figure 2-18 to meas- 
ure local ion currents collected by the chamber sidewall. We 
also documented the electron current collected by the propellant 
plenum and the ion flux to the cathode magnet and screen elec- 
trode under various operating conditions, and with these latter 
electrodes biased 27-V negative with respect to the cathode to 
ensure ion saturation. 

Figure 2-19 shows the distribution of electron and ion 
currents to the chamber boundaries and internal electrodes. The 
electron current is observed to be almost entirely collected by 
the two magnet rings located near the ion-extraction assembly. 
The ion current densities derived from the wall-probe data are 
only about 1/10 as large as the current density at the ion- 
extraction assembly, suggesting that ions preferentially drift 
to the downstream end of the chamber. 

When the wall probes were electrically connected to the 
cathode (through the current meters shown in Figure 2-10) their 
net current was positive (ion collection exceeded electron 
collection) and essentially independent of ion-production cost 
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Figure 2-17. Ring-cusp-thruster beam profiles. 
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Fiqure 2-18. Construction details of ion-flux probe. 
(Collector area is 0.74 cm 2 .) 




(cathode-emission current), as shown in Figure 2-20. The screen 
grid, on the other hand, was observed to change from a condition 
of net-ion current collection to one of net-electron current 
collection as the ion-production cost was increased. These 
measurements are also shown in Figure 2—20 , which indicates that 
the net-electron current tends to saturate at about ei - 180 
eV/ion. The transition from net-ion collection to net-electron 
collection occurs at about = 8.5 A, which corresponds to 
the "knee" of the performance curve (*>150 eV/ion for these 
conditions ) . 

6. Lanqmuir-Probe Measurements 

We documented the transition from net-ion to net-electron 
current collection by performing Langmuir-probe surveys under 
operating conditions that placed us on either side of the 
transition point. The diagnostic measurements were conducted 
using the probe-sweep and data-acquisition system shown in 
Figure 2-21 . The sweep circuit varies the probe voltage from 
-20 V to +80 V (with respect to cathode potential) over a period 
of about 1 sec. Isolation amplifiers are used to isolate the 
probe circuit from the data-acquisition system, which consists 
of a multiplexed, 16-channel analog-to-digital converter 
interfaced with an S-100 Microsystems computer. The digital 
results are stored on floppy disks, with the information later 
used by a VAX 11/780 computer that performs the numerical 
analysis 2-1 * leading to the plasma properties. 

interpretation of the probe characteristics showed that the 
electron-energy-distribution function is predominantly 
Maxwellian beyond a central "core" which extends downstream of 
the cathode. Within the core, the distribution function is the 
sum of a monoenergetic (or primary) plus a Maxwellian contribu- 
tion (see Appendix A). The spatial variation of the plasma 
properties shows that the two-group description fits the 
measurements quite well throughout the ctoss-hatched region 
shown in Figure 2-22. In this sense, the spatial variation of 
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Figure 2-21 * Block diagram of probe-sweep and data- 
aquisition system used for conducting 
Langmuir- probe measurements. 
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Figure 2-22* Region of the ring-cusp discharge chamber 
where the two-group-p.i asma assumption is 
valid. 


54 



of the electron-energy-distribution function within the ring- 
cuso discharge chamber is similar to those observed in other 
configurations (such as the divergent-field thruster); the 
primary electrons exist predominately in the central region 
defined by the location of the cathode and the "critical 
magnetic-field line." Beyond this boundary the electron energy 

and ionization rate drop rapidly. 2 ” 5 

An indication of the location of the gradients in plasma 
properties is evident in Figure 2—22 in the form of the density 
of the "dots" , which denote the locations at which we performed 
our final probe measurements. An initial survey was conducted 
in which we determined the minimum number of probe locations 
required to characterize the plasma by performing measurements 
at each of the locations indicated by the mesh points of 
Figure 2-22. The closeness of the "dots" representing the 
probe— measurement locations is an indicator of the regions where 
steep gradients in plasma properties could be detected by 
comparing the probe characteristics. Figure 2-22 also shows the 
5-eV Maxwellian-electron isotherm, which is seen to roughly 
coincide with the boundary of the "core" region of the plasma, 
separating the uniform and non-uniform regions. 

Viewing the ring-cusp discharge chamber as others (such as 
the divergent field and single-cusp configurations) have been in 
the past allows one to deduce that the primary ionization volume 
is bounded at one end of the chamber by the electron— emitting 
cathode, and at the other end by the cathode-potential screen 
grid. The remaining boundary is established by the surface of 
revolution defined by the magnetic-field lines that pass near 
the cathode orifice and terminate on the magnet rings. This 
description leads to a predicted ion-production volume that has 
the approximate geometry shown in Figure 2-23. 

The electron- temperature distribution is shown in Fig- 
ure 2-24. The isotherms have nearly the same shape as the 
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magnetic flux lines 2 " 6 ; their divergence into the volume 
downstream of the cathode results in relatively high (« 3 eV) 
temperature in the region near the screen electrode. The 
electron-density contours have the same general shape as the iso- 
therms. The ionization collision frequency (defined as 
v = n Q + , where n Q is the electron density and Q q + is the 
temperature— dependent ionization reaction rate) is presented in 
contour form in Figure 2-25 (the ionization collision frequencies 
have been normalized to the maximum value, v max ), again showing 
the same ;eneral shape 2 " 5 as the magnetic field lines. The 
strona gradients in plasma properties are especially evident in 
this form, which shows that most of the ion-production occurs in 
the downstream region of the chamber, and that the bulk of the 
volume in the region upstream of the cathode contributes little 
or nothing to the ionization process. Plasma-potential contours 
are presented in Figure 2-26, showing the presence of a favorable 
potential gradient in the region upstream of the screen 
electrode, where most of the ionization occurs. 

For purposes of comparison. Figure 2-27 presents collision- 
frequency contours for a 30-cm-diameter state-of-the-art J-series- 
equivalent thruster that were obtained under NASA contract 
NAS 3-21040. The contours diverge in the downstream direction, 
consistent with the divergence of the magnetic field. The 
spatial extent of the ion-production region in the upstream end 
of the chamber is noticeably greater than it is in the ring-cusp 
chamber (Figure 2-25) due to "spreading" of the source of 
electrons by the baffle and polepiece. The Maxwell ian-electron 
isotherms have a shape similar to the collision-frequency 
contours, but the temperature near the center of the screen grid 
in the J-series-equi valent thruster is substantially lower than 
it is in the ring-cusp configuration. 


*To simplify the calculations, we neglected the contribution 
due to the primary-electrons and used the plasma properties 
obtained under the assumption of a strictly Maxwellian 
distribution function. 
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Fiqure 2 - 26 . Plasma-potential contours in the 
mercury ring-cusp thruster. 
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The results presented in Figures 2-24 through 2-26 corre- 
spond to operation near the "knee" of the performance curve, 
which in the discussion of boundary currents was indicated to be 
characterized by net-electron current collection at the cathode- 
potential screen electrode. Plasma-property measurements at 
lower eV/ion (where net-ion current is collected) have the same 
general characteristics as those obtained at higher eV/ion. How- 
ever, the electron temperature is substantially lower when the 
thruster is operated under the less-efficient discharge condi- 
tions. This effect is evident by comparing Figure 2-24 with 
similar results presented in Figure 2-28. The overall tempera- 
ture is higher and the isotherms are much more divergent under 
the high-eV/ion operating conditions. The temperature near the 
center of the screen grid is > 8 eV, allowing electrons in the 
high-energy tail of the distribution function to reach the 
cathode-potential surface. Probes located downstream of the 
accelerator electrode (no ion— beam extraction) in an inert— gas 
thruster 2- 1 show that the collection of net— electron current 
occurs only in the center of the ion— extraction assembly, 
supporting the argument presented above. 

7 . Preliminary Model of the Ring— Cusp Discharge Chamber 

Using the results presented above, we can formulate the 
following (partial) description of the plasma processes occurring 
within the ring-cusp discharge chamber: 

• The bulk of the plasma is concentrated in the 
downstream end of the discharge chamber by the 
divergent-cusp magnetic field. 

• Electrons emitted by the cathode are "guided" into the 
ion-production volume by the magnetic-field lines 
passing through the region near ’.he cathode orifice. 
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» Electrons are repelled from the chamber (anode) 

boundaries by the large values of the integral JBdA, and 
are predominantly collected by the cusp regions of the 
magnets located in the downstream end of the discharge 
chamber. 

9 Ions are preferentially lost to the downstream 
boundaries 

In conventional divergent-field, single-cusp, and multipole 
discharge chambers, the magnetic fields are on the order of 100 G 
and, therefore, have little effect on ion trajectories. However, 
with magnetic fields on the order of 1,000 G typical of the ring- 
cusp thruster, the ion motion can be altered. Kaufman ^ ® has 
formulated the magnetic— f ield integral necessary to repel an elec- 
tron of energy, c p (eV), from a magnetically shielded anode as 

/ BdA > 6.74 , G-cm , (2-2) 

where /Bd* is the integral of the magnetic-field component per- 
pendicular to the path £. Using a similar approach, the magnetic- 
field integral required to repel a mercury ion travelling with the 
3ohm velocity* in a direction normal to a magnetically shielded 
surface can be formulated as 

/Bd* > 1442 v/T^r G-cm , (2-3) 

where T m is the electron temperature in eV. The results 
obtained from Equations (2-2) and (2-3) are plotted in I ig- 
ure 2-29 , showing that for typical plasma conditions of ? p = 30 
eV and T m = 5 eV, the corresponding magnetic-field integrals are 
about 40 and 3,000 G-cm for electrons and mercury ions* 
respectively. For argon ions, the integral drops to about 1,400 
G— cm because of the reduction in atomic mass. 

*This approach is conservative, since there is evidence 2 
that ions approach magnetically shielded anode-potential surfaces 
at a rate substantially less than that corresponding to Bohm 
diffusion. 
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Using the magnetic-field measurements described above, we 
calculated values of the integral / Bd£ along various paths to the 
chamber walls. Our results are shown in Figure 2-30, indicating 
that values of the integral on the order of 1,000 G-cm are typi- 
cal of this design. We believe, therefore, that ion trajectories 
can be influenced by the strong magnetic fields established by 
the ring cusp confinement scheme. The magnetic deflection scales 
as the square root of atomic mass, so that lighter ions such as 
argon would be influenced to a greater degree. This should lead 
to better performance with lighter gases, a prediction that is 
borne out by the experimental results 2 " 7 obtained with a ring- 
cusp geometry similar to the one investigated here, and which 

performed somewhat better (25% lower eV/ion, using argon as the 
propellant) . 

B. SINGLE-CUSP DISCHARGE CHAMBER 

In addition to the high-performance ring-cusp discharge 
chamber, we also investigated two alternative designs (the 
single-cusp and multipole) that employ weaker, boundary-type 
magnetic fields to produce a relatively field-free plasma volume, 
which should permit full utilization of the lifetime and 
perveance capabilities of the ion-extraction assembly. The 
single cusp— magnetic— field configuration is conceptually similar 
to that of the ring-cusp arrangement, but with fewer cusps and 
weaker fields. The multipole-magnetic-field configuration is 
similar to the single-cusp arrangement, but with an increase in 
the number of cusps. In both the ring-cusp and single-cusp 
configurations, the electrons emitted by the cathode have access 
to the plasma volume bordering the ion-extraction assembly by 
transport along the diverging magnetic field lines, much like 
the conventional divergent-field design. By contrast, the elec- 
trons in the multipole-magnetic-field arrangement are emitted 
into a weak-f ield region, with the strong-field region near the 
chamber sidewall and endwall used to repel the energetic 
electrons. 
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Figure 2-30. Measured values of the integral / Bd i 
in the ring-cusp discharge chamber. 
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A sketch of the laboratory-model single-cusp thruster is 
shown in Figure 2-31. The discharge-chamber configuration is 
arrived at by incorporating a relatively minor structural 
modification in the state-of-the-art J-series thruster (the 
addition of an upstream anode and pole piece) to produce the 
desired magnetic-field geometry. Axial and radial electromagnets 
are employed in the laboratory-model configuration, permitting 
rapid variation of the magnetic-field strengths. Separate power 
supplies (and a switching arrangement) enable the magnetic-field 
geometry (in the operating thruster) to be changed from the 
conventional divergent-field geometry to the advanced single-cusp 
configuration. This feature allowed the effects of a radical 
change in plasma-confinement technique to be evaluated in real- 
time, while maintaining essentially constant thruster-operating 
conditions. 

The discharge chamber was equipped with the movable cathode 
assembly shown in Figure 2-32. We designed and fabricated this 
highly functional apparatus, and we utilized it to optimize the 
locations of the cathode and associated components in each of the 
discharge-chamber configurations that we investigated (the ring- 
cusp, single-cusp, and multipole). The movable-cathode 
arrangement of Figure 2-32 enabled the axial position of the 
cathode pole piece, magnet pole piece, baffle, and cathode-keeper 
assembly to be independently varied during thruster op ration. 
This is accomplished using a single drive shaft to engage one of 
the four drive gears that transfer motion to the corresponding 
cathode component. The appropriate gear is selected using a 
"gear-shift" arrangement, and the location of each component is 
tracked and displayed using a potentiometer and a digital 
readout. A solenoid positioned around the outside of the cathode 
pole piece was used to produce the magnetic field within the 
baffle-aperture region; the open area of this region could be 
varied by replacing the (movable) baffle plate. The pole piece, 
baffle, and baffle supports were fabricated using low-sputter- 
yield niobium. - 
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igure 2-31 . schematic of Hughes 30-cm 
single-cusp thruster 
(shown with high-emission- 
current cathode described 
in Section 4). 






Figure 2-33 presents performance measurements obtained with 
this thruster (operated as both a conventional divergent and 
advanced single-cusp configuration), showing that good 
performance (comparable to that of the highly optimized J-series 
thruster) was obtained with both magnetic-field arrangements. A 
comparison of the beam-current-density profiles obtained with the 
single-cusp and divergent magnetic-field arrangements is 
presented as Figure 2-34. The single-cusp profile is much more 
uniform, as anticipated. The reduction in maximum current 
density results in a 45% increase in the beam— flatness parameter, 
which implies a corresponding increase in the lifetime of the 
screen electrode. 2 * The calculated flatness parameter of the 
single-cusp profile of Figure 2-34 is F = 0.85, based on the 
active area of the electrodes. This high degree of uniformity is 
comparable to the results preseated below for the multipole 
thruster. 

C. MULTIPOLE DISCHARGE CHAMBER 

The design details of the 30-cm-diameter multipole thruster 
are shown in Figure 2—35. The sidewall of the discharge chamber 
uses 2.7-cm-wide sections consisting of the magnets, pole pieces, 
anode, and insulators. These sections are "stacked" together to 
produce the desired overall chamber length. The criteria we 
followed in designing the chamber was to provide an integral of 
magnetic-field strength (along a path perpendicular to the 
surface of the anode) of sufficient magnitude to deflect 
energetic electrons away from the anode. Following the work of 
Kaufman, 2 ” 8 a value of the integral / Bd* equal to 60 Gauss-cm is 
sufficient to "repel" a 60-eV primary electron. The multipole 
discharge chamber was initially set up to provide about 1.5 
times this value of /Bd £# and through the use of overlapping, 
variable-length anodes# we could increase or decrease this value 
to optimize the discharge-chamber performance. During prelim- 
inary tests we (iteratively) re-positioned the anodes (relative 
to the pole pieces) to provide nearly equal anode current 
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Figure 2-33. Performance measurements obtained with 
the single-cusp thruster operated with 
two magnetic-field geometries 
(J b = 2 A, V D = 32 V). 



DIVERGENT 



^ Figure 2-34* Beam profiles obtained with the single-cusp 

thruster operated with two magnetic-field 
t geometries (J b = 1.3 A, V Q = 32 V). 






densities, while at the same time satisfying the requirement of 
the magnetic-field integral. An example of the results of such 
an iteration is presented in Table 2-3. 

We evaluated the performance of the 30-cm-diameter multipole 
discharge chamber over a range of length-to-diameter ratios by 
removing sections of the cylindrical sidewall, as shown in Figure 
2-35. The performance variation is shown in Figure 2-36, and the 
corresponding beam-current-density profiles are presented in 
Figure 2-37. The performance data show the anticipated reduction 
in maximum propellant utilization as the chamber length (neutral 
residence time) is reduced. Little reduction in baseline eV/ion 
was realized as a result of shortening the chamber, suggesting 
that the sidewall losses were not as significant as the endwall 
and screen-electrode web losses. 

We measured ion losses to the chamber anodes (which comprise 
most of the interior area) by biasing them (individually) 20-V 
negative with respect to the cathode. The results are presented 
in Table 2-4 for the discharge-chamber configuration of Figure 2- 
35(c), which has a length-to-diameter ratio of L/Dfc = 0.18. 

The results are presented as current densities, where we used the 
width of the anode to define the ion-loss area. Table 2-4 also 
presents the discharge current densities, where the area between 
adjacent pole pieces was used to define the electron-collection 
area. The average ion flux to the screen electrode is given by 


where Afc is the beam area (641 cm 2 ), and <f> s is the ion 
transmission of the screen grid. Using <j> 3 = 0.674 (geometrical 
transmission) gives j s = 3*3 mA/cm 2 for a beam current of 
Jb * 1.5 A. A comparison of this value with the measured 
results of Table 2-4 indicates a significantly higher ion flux 
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Table 2 - 3 . Normalized Anode Current Densities in the 
Multipole Thruster (Anode Numbering 
Scheme is Shown in Figure 2 - 35 ) 




ANODE NO. 


Before 

Repositioning 

Anodes 


After 

Repositioning 

Anodes 







DISCHARGE PROPELLANT-UTILIZATION EFFICIENCY, r? m , % 


Figure 


36. Performance obtained with the multipole 
thruster configurations shown in 
Figure 2-35. 














Table 2-4. Discharge and Ion Current Densities in the 
Multipole Thruster Showing the Effect of 
Anode Position (L/D 5 = 0.18) 


Current Density, mA/cm 2 

Anode Before Recessing Anode #3 After Recessing Anode #3 


Electron 

Ion 

Electron 

Ion 

6.78 

1.65 

7.37 

1.66 

5.94 

2.21 

6.96 

2.24 

15.84 

3.04 

6.94 

2.00 

8.92 

1.98 

10.49 

2.07 

21 .84 

2,79 

24.84 


14.4 

1 .88 

15.05 

1.99 







density at the screen electrode.* Comparing the ion current 
arriving at the screen electrode (Jb/'t's) with the sum of the 
measured ion currents indicates that about one-half the ions 
produced in the discharge arrive at the screen electrode even 
though it comprises only about one-third of the total bounding 
area. 

Both the electron and ion-current-density measurements 
suggest that anodes #3 and #5 act as a plasma sink. To verify 
this, we recessed anode #3 to see if any performance improvement 
could be realized. Table 2-4 presents the corresponding current- 
density measurements, showing a substantial reduction in both the 
electron and ion losses for anode #3. Inasmuch as the baseline 
discharge losses for this configuration were about = 225 
eV/ion, we concluded that recessing anode #5 tc achieve more 
uniform electron and ion losses would probably not improve the 
performance in any significant manner. Considering the high 
performance demonstrated by the ring-cusp discharge chamber, no 
attempts to optimize the performance of the multipole 
configuration were made beyond those described above. 

Although the performance of the multipole discharge chamber 
was somewhat below our calculated predictions, some valuable 
information related to preferential ion drift and sensitivity to 
chamber length was obtained. In addition, it was shown that the 
conductivity in the baffle-aperture region can be controlled 
using a magnetic-baffle arrangement that does not require an iron 
circuit. In the cathode pole piece/baffle design of Figure 2-32, 
the magnetic field in the vicinity of the baffle aperture is 
produced by the fringe field of the external coil. Figure 2-38 
presents measurements showing the effects of magnet-coil current, 
J mb' on cathod ®~vaporizer temperature, T cv . These results, and 


*More recent results obtained by Kaufman 2-9 and Wilbur 2-10 
have shown that for a variety of ion-chamber configurations, the 
ions preferentially drift toward the cathode-potential surfaces 
(as compared to magnetically shielded anode-potential 
surfaces) . 




the calibration of Figure 2-14(b), allow one to infer that the 
fringe-field baffle arrangement can be used to vary the cathode 
flow rate over a wide range (to permit discharge-voltage control 
or to accommodate throttling). 

Perhaps the most outstanding performance characteristic of 
the multipole discharge chamber is the attainment of highly 
uniform plasma conditions within the relatively magnetic-field- 
free volume of the plasma. This uniformity is manifest in an 
extremely uniform beam-current-density profile. Figure 2-39 
presents the beam-current-density variation for the multipole 
thruster operated at a beam current of J5 = 4 A. The flatness 
parameter for this profile based on the active (apertured) area 
of the electrodes is F = 0.95, which we believe to be the highest 
value that has been reported to date. 
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Figure 2-39. 


CHAMBER RADIUS (ARBITRARY UNITS) 


Beam-current-density profile of the multipol 
thruster operated at J b = 4 A ( L/D B * 0.36). 



SECTION 3 

ION- EXTRACTION- ASSEMBLY EVALUATION 


Hiqh oerveance ion-extraction assemblies require an 
electrode-mountinq atranqement that is dimensionally stable over 
a wide ranqe of input (discharge-chamber) power. Thermo- 
mechanical analyst l has shown that the "cold" interelectrode 
soacinq of the state-of-the-art J-series ion-extraction assembly 
shown in Fiqure 3-1 (a) decreases during ©Deration as a result of 
thermal expansion of the electrodes and mounting rings. The 
analysis also showed that flexible "leaf soring" supports permit 
the screen- electrode and mounting-ring assembly to freely expand, 
resulting in less interaction with the optics mounting ring and 
an increase in the dish depth of this electrode. The use of 
rigid posts for supporting the accelerator-electrode mounting 
ring, on the other hand, results in a bending moment that deforms 
this mounting ring and inhibits downstream displacement of the 
electrode. This same study showed that the flexible accel- 
support-ring mount of Figure 3-1 (b) allows the electrode and 
support ring to freely expand, increasing the dish depth and 
resulting in reduced sensitivity of the interelectrode spacing to 


variations in discharge power. 

Extracting high beam currents at low specific impulse can be 
accomplished by the addition of a decelerator electrode to enable 
the ratio of net-to-total accelerating voltage to be reduced to 
as low as R * 0.2 without experiencing excessive beamlet diver- 
gence. The three-grid ion-extraction assembly of Figure 3-1(0 
employs flexible mounting-ring supports for all three electrodes 
and was developed for the high-discharge-power and low-hetm-power 
operation required of the advanced ion-cnruster designs. Both 
advanced-technology ion-extraction assemblies of Figure 3-1 were 
investigated under this program. The electrode-aperture 
geometries are indicated in Figure 3-2. 
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(a) TWO-GRID OPTICS S/N 917 



(b) THREE-GRID OPTICS S/N 914 


Figure 3-2. Electrode geometries of the ion- 
extractions assemblies that were 
evaluated. (Dimensions in inches.) 



A. 


TWO-GRID ASSEMBLY 


The oerveance of an ion-extraction assembly should be inde- 
cendent of thruster operating conditions as long as the ixiter- 
electrode soacing remains constant and the current density 
oroflle does not change. However, measurements conducted using 
the ion-extraction assembly of Figure 3-1 (b) and a conventional 
J-series-eauivalent thruster (S/N 301J) are presented in Fig- 
ure 3-3, showing that that the perveance increases with discharge 
voltage. The measurements were oerformed under conditions of 
constant discharge oower, and, therefore, the perveance variation 
cannot be attributed to power-induced changes in interelectrode 
soacing. Instead, it is attributed to the increase in multiple 
ionization that accompanies the improvement in oropellant- 
utilization efficiency that is characteristic of operation at 
higher discharge-voltage and beam-current levels. The analysis 
presented in Appendix B shows that the measured perveance should 
be corrected for the presence of multiply charged ions since the 
theoretical perveance scales as the square root of the charge-to- 
mass ratio of the ion. The relationship between measured and 
corrected perveance is shown to be 

^corrected ^measured ’ (3-1) 

where P is the perveance, and a is the thrust factor that cor- 
rects measured thrust for the presence of multiply charged ions. 
Usinq this relationship to correct the measured perveance we 
generated the lower curve of Figure 3-3 which shows that the 
corrected perveance is invariant with discharge voltage.* These 
results indicate that meaningful perveance comparisons should be 
made only after correcting the measured results for the effects 
of multiply charged ions. 

*The discharge-volt age- independent results of Figure 3-2 were 
obtained using the average value of the correction factor a. 

The centerline or peak value "over corrects" the perveance 
results, consistent with the results presented below, showing 
that the location of the perveance limit is off-axis. 
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PERVEANCE 




1- ~gure 3-3. Variation of perveance with 

discharge voltage for thruster 
S/N 301J equipped with the 
advanced-technology ion- 
extraction assembly S/N 917 
shown in Figure 3-1 (a). 
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Figure 3-4 presents the variation of accelerator drain 
current with total accelerating voltaae and beam current, and the 
oerveance lines (corrected for doubly charged ions) derived from 
these data. The change in oerveance (30^) over the beam-current 
ranqe 0.75 A < J^ < 2.0 A is less than half that observed with 
the state-of-the-art ion-extraction assembly, indicating that the 
interelectrode spacing of the advanced ion-extraction design 
could be less sensitive to variations in discharge power. As a 
result of maintaining nearly constant interelectrode spacing, the 
oerveance of the advanced ion-extraction assemblies should be 
less sensitive to discharge power, permitting stable operation 
over a wider range of input power. Starting with identical 
"cold'' interelectrode spacings, the J-series ion-extraction 
assembly can be expected to have higher oerveance (as compared to 
the advanced ion-extraction assembly) as the power is increased 
and the interelectrode spacing becomes smaller. This expectation 
was verified by comparing the minimum-total-voltage requirements 
of the advanced ion-extraction design with those representative 
of the J-series thruster. Figure 3-5 presents the comparison, 
which shows that the advanced-technology ion-extraction assembly 
exhibits comparable or hiqher perveance for beam currents up to 
Jb = 1.3 A. At higher current levels, the perveance of the 
J-series ion-extraction assemblies is greater; this could be an 
indication of a larger reduction in their interelectrode 
spacings. The lower perveance of the f lexible-accel-support 
design should, however, be offset by an increase in its maximum 
current-extraction capability which results from the improved 
dimensional stability. 

























Figure 3-6 presents perveance measurements that were 
performed with the multipole thruster described above, and using 
the rigid and flexible mounting arrangements of Figures 3-1 (a) 
and 3-1 (b) to support the same pair of electrodes (with the same 
"cold" interelectrode spacing). The results of Figure 3-6 show 
that the perveance of the flex-mount ion-extraction assembly is 
greater than that corresponding to the rigid-mount arrangement. 
This result is contrary to the stress-analysis results which 
predicted that the interelectrode spacing of the flexible-mount 
arrangement would remain nearly constant, while that of the 
rigid-mount arrangement would decrease with increasing power 
level. The results of Figure 3-5 were obtained with a J-series- 
equivalent thruster (peaked plasma-density profile), while those 
of Figure 3-6 were obtained using the multipole thruster (uniform 
plasma-density profile). Therefore, we speculate that the 
thermally dependent perveance characteristics are also a function 
of the local plasma-density (and power) distribution. 

The results of our measurements suggest that the location of 
the perveance limit is not necessarily at the location of highest 
current density (center of the electrodes). The bases for this 
argument are the results discussed above, plus two additional 
measurements in which we showed that ( 1 ) the perveance of an ion- 
extraction assembly that was operated with uniform and nonuniform 
beam profiles was higher for the more peaked profile, and (2) the 
perveance of an ion-extraction assembly was independent of the 
aperture displacement at the boundary of the electrodes. In the 
first measurement, the single-cusp thruster was operated at 
constant discharge power and beam current, but with the magnetic 
field arranged to produce either the uniform or peaked beam- 
current-density profiles of Figure 2-34. Contrary to our 
expectations, the flatter profile exhibited lower perveance, 
suggesting that the perveance limit occurs off-axis,* i.e., beyond 


ACCELERATOR CURRENT, J 



Figure 3-6. Performance of ion-extraction assembly 
S/N 917 operated on the multipole 
thruster. The electrode mounts were 
changed from the rigid ones shown in 
Figure 3-1 (a) to the flexible ones 
shown in Figure 3-1 (b ) , while 
retaining the same electrodes. 
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the radial location defined by r = 10 cm, where the current den- 
sity of the uniform profile exceeds that of the nonuniform dis- 
tribution. In the second measurement, we performed permeance 
measurements using an electrode set that had substantial aperture 
offset in the region defined by r > 12 cm. The (unintentional) 
aperture misalignment was large enough to cause elongation of the 
accelerator apertures due to direct impingement of high-energy 
ions. Accelerator-current measurements obtained with this elec- 
trode set are presented in Figure 3-7. The measurements were 
performed about two-hours apart and indicate that the baseline 
accel current was decreasing at the rate of about 1 mA/hr due to 
elonqation of the apertures located at the periphery of the 
accelerator grid. If, however, the baseline accelerator current 
is subtracted from both sets of measurements, the resultant data 
fall on the same curve (Figure 3-7); therefore, the perveance is 
unchanged (independent of edge-aoerture geometry, which did 
change). Since the perveance limit corresponds to a condition of 
excessive accelerator current, we concluded that the radial loca- 
tion of the perveance limitation had to be less than that corres- 
ponding to the edge apertures which were already operating under 
conditions of excessive impingement current.. 

Other measurements 3-3 obtained using a single-aperture 
electrode arrangement with aperture displacement and. current 
representative of the electrode compensation and plasma- density 
profile of the J-series thruster showed that the location of 
the perveance limit is in the region 9 cm < r < 12 cm. These 
findings suggest that an increase in the perveance of the J- 
series ion-extraction assembly could be achieved by eliminating 
the aperture compensation, or by increasing the accelerator- 
aperture diameter in the region 9 cm < r < 12 cm. 

B. THREE-GRID ASSEMBLY 

We conducted both an experimental and analytical investiga- 
tion of three-grid ion-extraction assemblies. The experimental 
study was carried out using a government-furnished ion-extraction 
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Figure 3-7 . Variation of accelerator current with total 
accelerating voltage for electrodes that 
had poorly aligned apertures in the region 
near the periphery. 




assembly that employs the design features shown in Fig- 
ures 3-1 (c) and 3-2 (b). The analytical study was carried out 
using an adaptation of a three-electrode beam-steering model that 
was develoraed by Conrad . 3 " 2 

1 . Experimental Investigation 

Extracting beam currents as high as Jfc = 6 A at low 
specific impulse (low beam voltage V 5 ) requires a decelerator 
electrode to position the neutralization plane close to the accel 
electrode, enabling the net-to-total accelerating-voltage ratio 
to be 1 educed to as low as R = 0.2 without experiencing excessive 
beamlet divergence. The three-grid ion-extraction arrangement of 
Figure 3-1 (c) was designed for the high-beam-current and low- 
specific- impulse operation required of the extended-performance 
thruster; flexible mounting-ring supports are used for all three 
electrodes. Extrapolating the results of the thermomechanical 
analysis 3 " 1 described above indicates that this arrangement will 
allow the electrode-supporting rings to expand radially as they 
are heated by the discharge, maintaining nearly constant 
interelectrcde spacing. 

The three-grid ion-extraction assembly was operated with the 
ring— cusp discharge chamber for beam currents up to j)^ = 5.3 a. 
Measurements showing the variation of accel and decel currents 
with total accelerating voltage are presented in Figures 3-8 and 
3-9 for beam currents in the range 2 A < J b < 5 A. Figures 3- 
10 and 3-11 present the measured variations of accel- and deeei- 
electrode currents with net-to-total accelerating-voltage ratio, 

R. These results* demonstrated the ability to successfully 
operate cne 3-grid design of Figure 3-1 (c) over a net-to-total 
accelerating-voltage-ratio range of 0.2 < R < 0.8 for beam 
currents in the range 2 A < Jb < 5 A. 

*The aperture pattern in the ion-extraction assembly was under- 
compensated, resulting in excessive decelerator-electrode 
currents and noticeable elongation of the decel apertures near 
the- periphery of the electrodes. The design information 
presented in Section 3.B.2 should allow future designs to be 
properly compensated. 
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Figure 3-8. Variation of accelerator-electrode current 
with total accelerating voltage for the 
ring-cusp thruster equipped with the three- 
grid ion-extraction assembly S/N 914. 
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Figure 3-9. Variation of decelerator-electrode current 
with total accelerating voltage for the 
ring-cusp thruster equipped with the three 
grid ion-extraction assembly s/N 914. 






NET-TO-TOTAL ACCELERATING-VOLTAGE RATIO, R 


Figure 3-10. Variation of accelerator-electrode current 
with the ratio of net-to-total accelerating 
voltage for the ring-cusp thruster equipped 
with the three-grid ion-extraction assembly 
S/N 914. 
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E'igure 3-11. 


Variation of decelerator-electrode current 
with the ratio of net-to-total accelerating 
voltage for the ring-cusp thruster equipped 
with the three-grid ion-extraction assembly 
S/N 914. 





2 . Analytical Investigation 

Figure 3-12 presents the electrode compensation, 

A 

necessary to align the apertures along the radius of curvature, 
R c , as a function of the aperture separation, £ e s,c3 , defined as 




2A„ + t - + t, 
9 s / d a 


(3-2) 


The subscripts and superscripts s, a, and d refer to screen, 
accel , and decel, respectively, and t represents the electrode 
thickness and is the interelectrode spacing. In performing 
the calculations, we assumed that the accel-aperture pattern 
would not be compensated. In this case, the screen-electrode 
compensation, C^, is negative (screen apertures moved closer 
together), while the decel-aperture compensation, C^, is positive 
(decel apertures moved farther apart) and the amount of compensa- 
tion is determined by the following geometric relationship: 


r s,d 


50 





(3-3) 


where Rg is the radius of the active region of the electrodes. 

The results of Figure 3-12 are applicable to both 2-grid and 
3-grid ion-extraction assemblies; compensating by these amounts 
should produce beamlets leaving the accelerator system in a 
direction normal to its surface. For the J-series ion-extraction 
assembly the beamlet-normal compensation is calculated by 
Equations (3-2) and (3-3) to be C® = -0.183. 

The additional compensation necessary to vector the beamlets 
parallel to the thrust axis was calculated for 3-grid ion- 
extraction assemblies using a theoretical beam-steering model. 
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Figure 3-12. 


Electrode compensation required to align the 
apertures along the radius of curvature for 
the 30-cm-diameter thruster (29-cm beam 
diameter) . 





This theory relates the vectorinq anqle y to the displacements of 
the accel and decel aoertures throuqh the expression 


vs- — y _ ■ y (3-4) 

* 3y a 7 a 3T b Y b ' 

where y fi and y fa are the displacements of the accel and decel 
apertures, respectively, and ty/dy a , 3Y/3yb are the vectoring 
sensitivities, which are functions of the aperture geometry and 
the net-to-total acc<. lerating-voltage ratio, R. Figure 3-13 
shows the geometry and defines the nomenclature used in the 
model . 

In ion-propulsion applications it is desirable to use the 
same tooling for fabricating accel electrodes of both 2-grid and 
3-grid assemblies. Therefore, Equation 3-4 can be written for 
the case where the screen and accel apertures are displaced by y 

s 

and y^ (relative to the accel aperture), using the 
transformations 


y 


a 


(3-5) 


and 


+ • 


This gives 


Y 





(3-7) 
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Figure 3-13 


Geometry used in beam-steering model 



Figures 3-14 through 3-21 present the variation of vectoring 
sensitivities with net-to-total accelerating-voltage ratio R for 
various values of the aperture dimension, d # , and the parameters, 
and l 2 , given by 


and 

*2 = *d + fc d + “T '' (3-9) 

where the terms were defined in Figure 3-13. The calculated 
results indicate a much greater sensitivity of beam steering to 
accel-aperture displacement and show that the vectoring direction 
is opposite in sign for displacement of the decel apertures. 

Also, the effect of the decel electrode is diminished as R 
approaches unity (electric field between accel and decel elec- 
trodes approaches zero). Because of the difference in vectoring 
direction associated with displacement of the accel and decel 
apertures, maximum vectoring sensitivity is realized by dis- 
placing these apertures in opposite directions. Since the posi- 
tive displacement of the accel aperture corresponds to negative 
displacement of the screen aperture, maximum vectoring is 
achieved by displacing the screen and decel apertures in the same 
direction. Moving the decel aperture in the direction of beamlet 
vectoring should also reduce the possibility of direct 
interception on this electrode. 


0.150 


13111-7 


8 


d a ■ 0.045 in. 

a 


0.100 



-0.050 < I 1 1 I L 

0 0.2 0.4 0.6 0.8 1.0 


4 


2 


0 


2 


NET-TO-TOTAL ACCELERATING VOLTAGE RATIO, R 


Figure 3-14. Variation of vectoring sensitivities with R 
(d a « 45 mils). 
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Variation o£ vectoring sensitivities with R 
(d a = 60 mils) . 
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Figure 3-16. Variation of vectoring sensitivities with R 
(da « 75 mils ) . 
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The results of Figures 3-14 through 3-21 show that the 
effect of displacing the screen aperture is maximized by using 
large accel apertures* (large d^) and closely spaced electrodes 
(small and l 2 ) . The effect of displacing the decel apertures 
is also dependent on the interelectrode spacing and thickness, 
but, according to the model, it is independent of accel-aperture 
size. 

*The reduction in beamlet vectoring predicted for small accel 
apertures is consistent with the experimental observation that 
the J-series ion-extraction electrodes do not completely vector 
the beamlets parallel to the thrust axis, whereas the electrodes 
used in earlier designs (larger accel apertures) did. 


SECTION 4 


HIGH-EMISSION-CURRENT CATHODE STUDY 


In order to satisfy a cathode-emission-current requirement 
of J E > 30 A and a lifetime goal of L c > 15,000 hr requires 
that the state-of-the-art cathode design shown in Figure 4-1(a) 
e scaled for higher current operation. The basis that we used 
for scaling the cathode-insert design to a high-current version 
was to maintain conditions within the cathode that produce a low- 
voltage-arc discharge mode, and, as a result, efficient and lonc- 
ife operation. This discharge mode requires a certain minimum’ 
value of the product of pressure and characteristic dimension. *-l 
Therefore, the fundamental design consideration was to keep the 
pressure-distance product constant when scaling the dimensions of 
the J-series insert for use in the high-emission-current cathode. 
The characteristic dimension of the hollow cathode was assumed to 
be the inside diameter of the insert, dj.. 

The cathode-orifice diameter, d , should scale with emission 
current according to the relation 


where the exponent lies in the range 0.5 < n < 1. The lower 

value of n corresponds to constant-emission-current-density 
scaling, while the upper value of n is based on a correlation of 
experimental results. “-a Depending on the value of n. Equation 
(4-1) predicts that the cathode orifice diameter for J E * 40 A 
should be increased by a factor of 1.7 to 2.8 over the value used 
in the J-series cathode <J E . 12 A). Because increasing the 
orifice diameter would require an increase in cathode flow rate 
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Fxgure 4-1. Cross sections of the J-series and high- 
emission-current hollow cathodes. 




(to maintain constant-pressure conditions), we kept the inside 
diameter of the insert constant in the scaled design. 

The scaling considerations discussed above led to the high- 
emission-current cathode design shown in Figure 4-1 *b). The 
orifice plate is removable, which permitted us to evaluate 
different aperture diameters, d Q , using the same cathode/insert 
configuration. The insert was instrumented with thermocouples to 
enable us to assess the effects of various operating conditions 

on insert temperature. 

The high-emission-current cathode o£ Figure 4-1 (b) was eval- 
uated using the three aperture diameters listed in Figure 4-2. 

Me operated each cathode in a bell jar for a period of 500 hr at 
an emission current of J E - 40 A. A summary of operating condi- 
tions for these tests is presented in Table 4-1 . 

Figure 4-3 presents the variation of insert temperature with 
cathode emission current for both the high-emission-current and 
conventional cathode designs. The cathodes used a conventional 
(scaled) keeper and were operated in a diode arrangement using a 
water-cooled electrode to collect the emission current. A 
liquid-nitrogen trap between the bell jar and the mercury ^ 

diffusion pump maintained the bell-jar pressure at about 5 x 1CT 
Torr during the tests. These results show that the maximum 
insert temperature of the high-emission-current cathode (operated 
at J E “ 34 A) is only slightly higher than the corresponding 
insert temperature of the conventional cathode (operated at an 
emission current of J E - 12 A). A model"- 3 developed by Hughes 
under NASA Contract NAS 3-21040 relates cathode lifetime to 
insert temperature through the depletion of barium from the 
porous-tungsten matrix. Based on this model, the results of 
Figure 4-3 indicate that the lifetime of the insert used in the 
high-emission-ourrent cathode should be comparable to that used 

in conventional designs (= 30,000 hr). 

A scanning electron microscope (SEM) was used to examine the 

cathode orifice plates prior to and at the conclusion of each 
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Figure 4-3. Variation of cathode- insert 
temperature with emission 
current (d Q = 1.02 mm) . 




4 _i Time-Averaged Cathode Operating Conditions for 
Table 4-1. ^ e 2 ^ 00 g A . hr Tests of the High-Emission- 

Current Cathodes 


Cathode 

No. 

Orifice 
Diameter, d Q , 
mm 

Flow 

Rate, 

mA 

Keeper 

Voltage 

V 

Insert Temperature 

Anode 

Voltage, 

V 

Downstream, 
T lf # C 

Upstream * 
' c 

1 

0.76 

191 

6.3 1 

14.1 



2 

1.02 

238 

4.8 

14.0 

1178 

1096 

3 

1.27 

221 

5.0 

14.0 

1164 

1084 



O 
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500-hr test. Figure 4-4 presents photographs of the orifice 
olates, showinq the appearance of the apertures after the 20,000 
J\-h-_ tests. Substantial erosion of the 0 . 76-mm-diameter aoerture 
is apoarent in Fiqure 4-4 (a). The dimensions of the 1.0- and 
1 , 27-mm-diameter apertures appear to have changed very little, 
although the surface conditions do show the characteristic 
faceted appearance of a surface that has undergone ion 
souttering. Figure 4-5 presents a comparison of shadowgraphs 
taken of impressions of the downstream side of the orifice 
plates, indicating that the cylindrical throat region of the 
0.76-ram-diameter aperture was completely eroded away. The 1.0- 
rara-diameter aperture plate shown in Figure 4-5 (b) shows evidence 
of similar damage. The 1 . 27-mm-diameter orifice plate is shown 
in Figure 4-5 (c), with only slight evidence of aperture erosion 
apparent. The dimensions of this aperture were measured after 
the test, and they were within the specified fabrication 
tolerances. 
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(a) d Q - 0.76 mm (78° VIEW, 137 x MAGNIFICATION) 



(b) d Q = 1.02 mm (60° VIEW, 71 x MAG N IF 'CATION) 



(c) d Q - 1.27 mm i60° VIEW, 65 x MAGNIFICATION) 

Figure 4-4. 

Photographs showing the upstream sides 
of the cathode orifice plates and the 
appearance of the apertures after operating 
at an emission current of J_ = 40 A for 
500 hr. b 




e. d 0 «1.27mm 


Shadowgraphs of impressions made 
of the cathode orifice piates 
after operating at an emission 
Current of Jo ■ 40 A for 500 hr. 







SECTION 5 


ISOLATOR BREAKDOWN-VOLTAGE MEASUREMENTS 

To operate an extended-performance thruster at high beam 
current and high specific impulse requires beam voltages greater 
than the voltage-standoff capability of the J-series-thruster 
propellant-flow electrical isolator shown in Figure 5-1. 

However, scaling the isolator to a Agher voltage rating is 
readily accomplished by adding a sufficient number of isolating 
sections,* provided that the minimum breakdown voltage of 
individual sections is known. 

The breakdown voltage of adjacent pairs of isolator sections 
was measured using the instrumented seven-section propellant-flow 
electrical isolator shown in Figure 5-2. The instrumented 
isolator was fabricated under NASA Contract NAS 3-21040, where it 
was used to measure the voltage variation along the length of the 
isolator using electrical leads that attach directly to the 
wire-mesh screens and exit through the isolator body. Figure 5-3 
is a schematic showing the setup used in conducting the 
breakdown-voltage measurements of individual sections of the 
isolator. The results are presented in Figure 5-4 where 
breakdown voltages are plotted as a function of propellant flow 
rate, ft. These results show that the breakdown voltage decreases 
with increasing flow rate, and, in general, increases with 
distance along the isolator. These trends are consistent with 
Paschen's Law and are a result of the pressure rise accompanying 
higher flow rates and the pressure drop along the length of the 
isolator. The results of Figure 5-4 indicate a conservative 
voltage rating of -340 V per section, suggesting a total voltage 
standoff of about 2,400 V. 


*The present design incorporates seven sections and has a 
conservative voltage-standoff rating of 1,500 V. 
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Fiqure 5-1 . P rope 11 ant- flow electrical isolator used in the 
J-series thruster. 
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Figure 5-3. Schematic of the test setup used 

for performing isolator breakdown- 
voltage measurements. 
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We also performed breakdown-voltage measurements for 
multiple isolator sections in order to verify that the total 
voltage-standoff is the product of the number of isolator 
sections and the breakdown-voltage per section. These tests were 
conducted at an isolator-base temperature of Tb ase = 320 °C, 
and at a vaporizer temperature corresponding to a flow rate in 
the range 450 mA < m < 500 mA. Figure 5-5 presents the measured 
breakdown voltages plotted as a function of the number of 
isolator sections, indicating a linear relationship between the 
two variables. All of our data tend to follow the same curve, 
indicating the absence of any "end effects." The slope of the 
line drawn through the data points is 330 V/section, which agrees 
with the section breakdown-voltage measurements of Figure 5-4 to 
within the precision of the measurements. 

The effect of varying the isolator-base temperature is 
illustrated in Figure 5-6, which presents the breakdown-voltage 
measurements for the sixth section of the seven-section isolator 
for base temperatures of 260°C and 320°C. Under low-flow 
conditions, a reduction in isolator-base temperature results in a 
substantial reduction in breakdown voltage. This effect is a 
consequence of Paschen's Law, which (in general form) relates the 
breakdown voltage, Vg, to the product of gas density, n, and 
interelectrode distance, d. The temperature dependence arises 
from the equation of state, resulting in n p/T for an ideal 
gas. This allows Paschen's Law to be formulated as 

v b -£[^) , (5-1) 

which reduces to the more familiar form, Vg * g(pd), for con- 
stant temperature, T. The discussion leading to Equation (5-1) 
suggests that the curves of Figure 5-6 should coincide if the 


■ ^ 


128 


BREAKDOWN VOLTAGE 


10810-10 


O ~ FEEDLINE GROUNDED 
□ ~ SECTIONS 1, 2 SHORTED 
A ~ SECTIONS 1, 2, 3 SHORTED 



NUMBER OF ISOLATOR SECTIONS 


Figure 5-5. Variation of breakdown 

voltage with the number of 
isolating sections. 
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Figure 5-6. 


Effect of isolator base 
temperature on the break- 
down voltage of the 
instrumented propellent- flow 
electrical isolator. 
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flow-rate data are divided by the absolute temperature of 
the qas, which can be assumed to be roughly the same as the 
isolator-case temperature. Figure 5-7 presents the results 
we obtained by replotting the data of Figure 5-6 versus the 

parameter *i/T bage and shows that the data points fit a sinc'e 
curve quite well. 



Figure 5-7. Variation of the isolator-section 
breakdown voltage with the 
parameter m/T haao . 



SECTION 6 


POWER-PROCESSOR SIMPLIFICATIONS 

Various approaches to reducing the number of power supplies 
required to operate the J-series thruster have been 
investigated, 6 " 1 ' 6-2 and it has been demonstrated 6 ” 1 that the 
number of individual power supplies can be reduced from twelve to 
five, using switching, combining, and power-absorbing techniques. 
Under this program, we investigated two approaches for 
simplifying the vaporizer power supplies of the J-series 
thruster, using efficient techniques that do not involve power 
absorption. 

The first concept uses a multiple-output power supply that 
provides for closed-loop control of the main, cathode, and 
neutralizer vaporizers, using the three switching regulators 
shown in Figure 6-1. This design requires three separate output 
filters (L x and C L ) which are the largest contributors to the 
volume and mass of the power supply. The remainder of the 
circuit consists of three high-speed switching transistors, Q 2 , 
three fast-recovery diodes, CR lf and three analog switching- 
regulator control circuits, U l . 

The second concept uses a time-shared or multiplexed power 
supply employing a switching regulator to sequentially provide 
power to each of the three vaporizers. A schematic of this 
design is shown in Figure 6-2. The multiplexed power supply 
eliminates many of the large and costly components of the 
multiple-output power supply by substituting low-level digital 
circuitry that is relatively small, lightweight, and inexpensive. 
This tradeoff is apparent in Table 6-1, which shows the parts 
required for each power-supply design. The multiplexed power- 
supply concept was selected as the most attractive candidate for 
further investigation because it offers a 70% reduction in parts 
count (and a corresponding increase in reliability). 
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We designed the multiplexed power supply of Figure 6-2 to 
orovide 5-msec current pulses to the cathode, main, and 
neutralizer vaoorizers at an operating frequency of 67 Hz. The 
cower supply was fabricated, and then evaluated using the J- 
series-equi valent thruster S/N 301J. The tests were performed 
using a two-inverter power processor that was supplied by NASA, 
with the multiplexed power supply used to replace its (three) 
individual vaporizer power supplies. 

Figure 6-3 presents a strip-chart recording showing the 
vaporizer-current waveforms of the multiplexed power supply. 

The vaporizer control loops are uncompensated, and we made no 
attempt to optimize their gain. Despite this, the power supply 
performed extremely well in these tests, and we were able to 
demonstrate the validity of the time-sharing concept. Fig- 
ures 6-4 through 6-6 show the response of the controlled 
variables (v D , J b , and V NK ) to step changes in the 
vaporizer currents that were obtained by momentarily reducing the 
output current to zero. The results demonstrate the stability of 
the unoptimized control system, and show recovery times on the 
order of < 40 sec (based on 1% deviation from steady-state 
values). Although we did not pursue the development further, we 
believe the system response and recovery time could be 
significantly improved through a modest optimization effort. 

The multiplexed vaporizer-power-supply concept represents a 
significant reduction in the mass and parts count of the three 
individual power supplies presently used in the power processor 
of the J-series thruster. This reduction is illustrated in 
Table 6-1 which compares the parts required for the J-series 
vaporizer supplies, along with the requirements for the Hughes 
multiple-output and multiplexed designs. The benefits are 
readily apparent by comparing the requirements for transformers, 
power transistors, and commutating diodes. The simplification is 
achieved at the expense of the low-level multiplexer logic 
circuitry. 
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Figure 6-3. 


Current waveforms of the multiplexed 
power supply operated at 67-Hz. 
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figure 6-6. “ 

output of the multiplexed power supply. 




SECTION 7 


THRUST- LOSS- MEASUREMENT ERROR ANALYSIS 

Accurate predictions of net thrust are required for 
planetary applications of ion— propulsion technology, since even 
small deviations could cause a significant loss in total impulse 
and pointing accuracy (when integrated along the flight path). 

The net thrust is obtained by correcting the thrust calculated 
from measured values of beam current, J^, and beam voltage, v fa , 
for losses resulting from beamlet divergence and multiply changed 
ions. These losses arise from (1) the reduced thrust-to-power 
ratio of multiply charged ions, (2) diverging ion trajectories 
caused by transverse electric fields within the acceleration 
region, and (3) vectoring caused by the non-planar shape of the 
electrodes and aperture displacement. The effects of multiply 
charged ions are accounted for in the factor a, while the effects 
of divergence and vectoring are accounted for in the factor F t . 
Both correction factors are less than unity (by definition), 
reducing the calculated thrust according to the relationship, 

F = aF t F c , ‘ (7-i) 

where F is the true thrust, and F c is the thrust calculated by 

F 0 - J b (2mV b /ej 1/2 , (7-2) 

where m is the mass of the ion, and e is the electronic charge. 
The accuracy of determining a, F t , and F c depends on both 
experimental and computational errors that are inherent in the 
measurement technique. This section describes an analytical and 


Pi?:iCL ? niN° r 


t:.i. 


r> 


141 


intentionally blank 


experimental investigation of the procedures that we use to 
determine the factors « and F t . a comprehensive examination of 
the measurement technique, assumptions, and equations used in the 
calculations leads to the identification of possible sources of 
error or uncertainty. An estimate of the overall accuracy of the 
thrust-correction factors is calculated using standard error- 

analysis techniques and estimates of the uncertainties in each of 
the variables. 


PROCEDURE FOR DETERMINING THE THRUST-LOSS FACTORS 

Obtaining the thrust-loss factors a and F t requires a 
measurement that can determine the singly and doubly charged ion- 
current-density vectors as a function of the polar coordinates r 
and $ defined in Figure 7-1 . Assuming symmetry about the thrust 
axis simplifies the measurement by eliminating the variable $, 
reducing the problem to that of determining the current density 
of each ion species as a function of the radial coordinate r. 

The measurement is performed using an articulating probe that can 
view a small region of the ion-accelerator electrode from differ- 
ent angles, 4 > p , and separate the total current leaving this 
region according to velocity (charge state). This gives the 
currents i n (r,* p ), where the index n is used to denote the charge 
state. Integrating these currents over the angle <j> p gives the 
magnitude of the current-density vector, T(r), and the angle * 
that this vector makes with respect to the thrust axis; 

Figure 7-2 illustrates the geometry. 

1 • Velocity Analyzer Probe 

The total current emanating from the active region of the 
accelerator electrode is separated into the current components, 
x n (r ' ^p*' using a series arrangement of a particle collimator 
and velocity filter. The collimator restricts the viewing area 
and transmits a highly collimated beam to the velocity filter. 


Figure 7-2. Illustration of the geometric vari- 
ables describing the current density 
at the accelerator grid. 




The transmission of the collimator is illustrated in Figure 7-3, 
which shows its response to a parallel beam inclined at an angle 
ij> 0 . The transmission is unity* when the collimator is aligned 
with the beam (<t> p = <t> 0 ) and rapidly drops to zero for angles 
about <|> 0 . The angle is the acceptance half-angle and is 
small enough ( «|> = 0.29°) to enable the probe response at any 
angle <|> 0 to be interpreted as the response due to particles 
that leave the viewing area and which follow straiqht-line 
trajectories inclined at an anqle (<j>p) with respect to the 
thruster axis. 

The velocity filter is composed of orthogonal electric and 
magnetic fields (E and 3, as shown in Figure 7-4) and transmits 
only those particles having velocity equal in magnitude to the 
ratio E/B. All other particles are "filtered out" by an 
imbalance in the electric and magnetic forces which deflects the 
particles up or down, depending upon whether the particle speed 
is greater than or less than the ratio E/B. The series arrange- 
ment of the collimator and velocity filter results in a probe 
output corresponding to collection of those particles having 
velocity vectors with magnitude E/B, and with direction in the 

range <|> p - < <t> < -*• -► 

A schematic of the velocity-analyzer or ExB probe that was 

developed 7-1 by Hughes for measuring i n (r,4> p ) is presented in 
Figure 7-5. The probe assembly consists of a collimator, 
velocity filter, drift tube, electron suppressor, and current 
collector. The collimator apertures have a diameter of 0.25 mm, 
resulting in an acceptance half-angle of ^ = 0.29°. With this 
qeometry, and with the probe positioned 38 cm downstream of the 
accelerator grid, the viewing area is 0.13 cm 2 , or the equivalent 
of about three accelerator aoertures in the J-series ion- 
extraction assembly. The separator provides orthogonal electric 

*The mathematical formulation and physical significance of the 
collimator transmission is described in Appendix D. 
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Figure 7-3. Angular response of the collimator 
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and magnetic fields E and £. The electric field is established 
by the potential applied to parallel plates, and the magnetic 
field is provided by a permanent magnet. Varying the plate 
potential changes the ratio, E/B, allowing ions of different 
velocities to traverse the separator undeflected, and to reach 
the collector. The collector aperture is biased 45 V negative 
with respect to the probe assembly in order to return secondary 
electrons to the collector. 

The dimensions of the collimator result in a small viewing 
area, but one which is large enough to provide a readily 
measurable current ( » 1 0“ 9 A) . The drift-tube length and 
collector-aperture diameter ensure collection of the entire 
collimated and undeflected beam of interest. Figure 7-6 presents 
collector current as a function of plate voltage and demonstrates 
the ability of the probe to resolve the peaks associated with 
singly and doubly charged ions; these measurements were obtained 
at an unusually high discharge voltage to exaggerate the Hg++ 
peak and to enable the Hg+++ peak to be resolved. The nearly 
rectangular shape of the peaks enables their height to be used as 
an indication of the total current of each species, simplifying 
the data acquisition and analysis. 

2. Test Facility 

A sketch of the probe setup used in the Hughes 9-ft-diameter 
vacuum chamber is presented in Figure 7-7. A precision stepping 
motor located on top of the vacuum chamber is used to move the 
probe vertically in or out of the thrust beam to vary the 
coordinate r. The pitch angle, $ p , is controlled using a 
precision stepping motor located inside the vacuum chamber, a 
potentiometer and digital readout provide a visual display of the 
probe pitch angle. The probe yaw angle can be adjusted to posi- 
tion the probe axis parallel to the thruster axis. The vacuum 
feedthrough is located off-center in the flange shown in 
Figure 7-7, enabling the lateral position of the probe to be 

adjusted (over a narrow range) so that its translational axis 
intersects the thruster axis. 
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Figure 7-6. Example of current output of the ExB 
probe as a function of deflection- 
plate voltage (with the probe 
inserted in the thruster beam) . — 
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Figure 7-7* ExB probe installation in Hughes 9-ft 
diameter vacuum chamber. 
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The only maintenance required of the orobe assembly is an 
occasional need to reolace the collimator aperture, which eventu- 
ally disintegrates as a result of sputtering by the ion beam. 

The absence of the collimating aperture is easily detected by the 
operator in the form of an increase in collector current and the 

inability to resolve the peaks corresponding to the individual 
ion species. 


3 . Thruster Alignment- 

The probe-to-thruster alignment is accomplished by aligning 
both the thruster and the probe with respect to the chamber axis. 
Thio approach minimizes setup time because the probe-to-tank 
alignment is required only after removal and reinstallation of 
the probe. The thruster-to-tank alignment is accomplished by 
aligning the thruster with respect to the axis of the vacuum- 
enclosure flange, as shown in Figure 7-8. The thruster support 
and vacuum-enclosure flange position the thruster with respect to 
the tank axis, and the angular alignment of the thruster axis is 
accomplished using the dial indicator to locate the edge of the 
accelerator grid at 90° intervals. In order to ensure that the 
annle between the thruster and flange axis is less than T, the 
difference between the readings obtained 180° aoart must not 
exceed 0.5 cm. if the difference is greater, the mounting fas- 
teners are loosened and the thruster is repositioned to obtain 
the necessary tolerance. 

4 • Data Acquisition 

The probe-positioning and data-acquisition system is semi- 
automated, requiring that the operator select the probe pitch 
angle, locate the beam edge, and determine the plate potentials 
corresponding to the various ion species.* The first step that 

of the 8- and 30-cm thrusters 

Therefore \K 9 ii!i5i% WUl !* 10n ° f triply char 9e<3 ions, 
rnerefore, the standard practice has been to measure emit, 

singly and doubly charged ion currents. The probe system does 

however, have the capability of measuring triply oh«gll ions 

beam h cu?ren?s. be ne9li * ible ln thrusters’ that^erate^t Mg"h' 
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the operator performs in conducting a beam scan is to set the 
probe pitch angle to <> p - 0, and then drive the probe down and 
into the beam until a collector current is detected. This 
locates the edge of the beam at the plane of the accelerator 
grid and provides a reference for subsequent movement of the 
probe into the beam, a total distance equal to one-half the 
diameter of the active region of the electrodes, if either the 
thruster or probe are removed from the test facility, then upon 
reinstallation the probe-to-thruster alignment is checked by" 
varying the probe yaw angle, a, until the probe current is at a 
maximum, with the probe-travel and thruster axes aligned, this 
step ensures that the viewing axis of the probe is coincident 
with the thruster axis. 

Next, the operator varies the output of the electric-field 
power supply to determine the voltages, v + and v ++ , corresponding 
to the peaks of the i + and i ++ current profiles (see Figure 7-4). 
These voltages are then input to the data-acquisition system via' 
potentiometers located on a control panel. Next, the probe is 
driven up and out of the beam, and the pitch angle, is set to 

the maximum negative value to be sampled (usually -15°). The 
probe is then driven down to locate the beam edge, which is 
always determined by moving the probe in the downward direction. 
This eliminates the possibility of a position error that could be 
caused by hysteresis in the positioning mechanism. At this 
point the operator places the system in the auto mode, causing 
the following steps to be executed 


• The probe is driven downward a distance equal to 
one-fourth the radius of the active region of the 
electrodes. 


• Any vibration of the probe caused by the directed 
motion is allowed to damp out. 
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® The currents i + and i ++ are sampled, printed, and 
punched on paper 'tape. 

The sequence described above is automatically repeated a total 
of four times, with the last measurement performed on the 
thruster axis. In the final step, the operator drives the probe 
up and out of the beam, increases the pitch angle by 5°, locates 
the beam edge, and resets the auto mode (which causes the steps 
outlined above to be repeated). 

In developing the computer program 7 " 1 used for analyzing the 
orobe data, one of the primary objectives was to minimize the 
number of incut data points, subject to the constraint that the 
accuracy of the reduced data should be at least comparable to 
that of the raw probe data (± 2% of full scale, based on the 
manufacturers specifications on the accuracy of the electro- 
meter). Minimizing the data-collection time ensures that the 
thruster operating conditions remain constant during the data 
scan, and relieves the operator of the time-consuming task of 
taking more data than is necessary. After the initial installa- 
tion of the probe, scans having a different number of data points 
were taken. It was determined 7 ” 1 that four equally spaced radial 
values and seven equally spaced angular values were near optimum. 
A larger number of radial data points produced essentially the 
same results, while a smaller number produced discontinuous- 
appearing curves. While six angular values for each radial value 
were usually satisfactory, seven were better for more divergent 
beams. The final technique chosen was to use four radial values 
and seven angular values, resulting in 56 data points per scan 
(28 each for i + and i ++ ). In practice, the operator may conduct 
the scan using as many as nine pitch angles, <j> p , so that the 
symmetry of the angular dispersion profiles can be checked, 
enabling the most symmetrical data to be used as input to the 
data-analysis routine. 
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3. DATA- ANALYSIS PROCEDURE 

The equations used to calculate the thrust-loss factors, a 
and F fc are derived in Apoendix C. The calculations are Performed 
using the magnitude, j(r), and direction, <{i(r), of the current- 
density vector at the ion— extraction accelerator grid. These 
variables are calculated from the ExB probe measurements, using 
the procedure derived and discussed in Appendix D. During the 
present investigation, we found discrepancies in the analysis 
procedure that had been employed to interpret the probe data, and 
we modified the procedure accordingly. Therefore, the descrip- 
tion presented here differs from that presented in the original 
reference 7 ” 1 describing the Hughes ExB probe. 

The angle <J> .is calculated at each radial coordinate using 
the expression. 



where i n (r, <f> p ) is the probe current corresponding to species 
n, <j> p is the probe angle, and <j> Q is the probe angle corresponding 
to the peak in the current profile. The magnitude of the current- 
density vector, j ( r ) , is calculated using the expression, 

4rrL 2 sec 3 <|> 0 f* max 0 

j n (r) = a~A_cos(y-4>) L i n (r ^p )cos V ec2( V 4> o ) l tan( V <t> o ) l d<|, p 


where A p and A Q are the probe- aperture area and probe- viewing 
area, respectively. The variable, L, is the distance from the 
probe to the accelerator electrode and is given by 

L = L© + R c (1 “ cos Y) f (7-5) 


A T ■ 9 
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where L is the distance from the probe to the electrode on the 
thruster axis, and R c is the radius of curvature of the 
electrodes (R c is positive for electrodes that are dished 
outward). The anqle y is qiven by 

Y = sin” 1 (r/R c ) * ^ 7 “ 


The total current density is obtained by adding the species 
current densities calculated from Equation 7-4; the total current 
of species n is calculated using the expression* 


J 


n 



n cos(Y 


*) 

* ' cos Y 


(7-7) 


where Rg is the radius of the active region of the accelerator 
electrode. The thrust-loss factors a and Ft ate calculated 
using the expressions 


o = 


, , £L _ 

J + * 2 J ++ 

J + + J ++ 


and 


(7-8) 


f (j + + 2^ COS * dA 

i_3 — 

11 (j t + f- n+) 


(7-9) 


F t = 


2» f Rg tj + + f- 1++) cos(T*) 

J o 

/ T 

J + + 2 


( 7 - 10 ) 


♦Equation (7-7) is obtained by combining Equations (C-1), (C-2), 
and (C-4 ) of Appendix C. 
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tin c*. 


.m, 


Integration of Equations (7-8) through (7-10) is performed using 
a second-order interpolation scheme, with the probe current 
extrapolated to zero to determine the integration limits <!> m ^ n and 

4>max • 

The ou tout of the data-analysis routine contains both 
graphical and numerical results. The plotted output consists of 
the angular dispersion profiles (a total of eight, corresponding 
to the two species and the four radial locations), and the radial 
current-density profiles (a total of three, corresponding to the 
two species and their sum) . The tabulated output consists of the 
values j,, j, , and $ (calculated at each radial location) and 

T TT 

the thrust-loss factors, a and Ff The tabulated results also 
contain the ratio of doubly to singly charged ion currents, and 
the propellant-utilization correction factor 3, defined as 


j. + j 

T 


(7-11) 


C. ERROR ANALYSIS 

The standard approach 7 " 2 used in evaluating the uncertainty 
U f , in a derived quantity f(x i ) is to calculate the sensitivity 
of f to each independent variable, xi , and then to estimate the 
uncertainty U Xi r in each of these variables. For equal odds on 
the uncertainties present, the overall uncertainty is given by 



where 3f/3xi is the sensitivity of the derived quantity, f, to 
the independent variable, xj.. 

The thrust equation is 

F = * (7-13) 
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and the uncertainty In the thrust can he written as 


u r 


/(£ °«f + (»; °» t ) + (»c u? c) 


(7-14) 


« u m assume that the uncertainty 
IO Simplify the discussion, we will assume 

in the calculated thrust is negligible. In this 
tion above simplifies to 

(7-15) 


-V p2 t °« + * 2 ° 2p t 


. . ,• lu 0 i-hrust can b0 GSfcitn&tsd 

, herpforA the overall uncertainty m the thrust ca 
rherefore, tn determined, 

if ^ uncertainties in the factors a and F t can be ae 

Because of the 0OmpU ^ ^“/Ipprl.irwM 0 ” determine the 

th e and 

to assess the sensitivity of the calculated re 

rrtamties in bbese wariables 

r„ ; almilar to that given by 

Equation q£ ^ f probe technique and ^equations 

used to analyse the probe output showed ^ ^ ^rdi- 
bles appearing in the expressions for a an t identifv a ll 

..... - v- 

the factors that could lead to sig the ExB measurement 

ties in these ^‘^he prtLry qua „tities listed in 

procedure r we . . could introduce uncertain- 

Tab le 7-1 as position variables that could ntr te 

ties into the variables, r and ♦ . Definitions of tn 

Uncertainties are provided in Figures 7-9 through 7-11. 


♦The expressions— r both 

fife" ^calculated results. 
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Radial coordinate. 

Angular coordinate. 

Electrode radius (active region). 

Inclination of plane of probe 
travel with respect to vertical 
diameter of thruster. 

Position of plane of probe travel 
with respect to vertical plane 
passing through the axis of the 
vacuum chamber. 

Probe yaw angle. 

Thruster pitch angle. 

Thruster yaw angle. 

Position of vertical diameter of 
thruster with respect to vertical 
plane passing through the axis of 
the vacuum chamber. 








r 

i 



i 

Figure 7-9. Definition of the coordinate uncertainties 

for: the variables r and <j>. - , 
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1 . Sensitivity Estimates 

We performed analyses to determine the sensitivity of the 
thrust-loss factors to the coordinates r and <j>p using the 
computer program described above to perform the numerical 
calculations. The baseline input data were taken from actual 
probe measurements obtained with a 30-cm J-series thruster 
operated at = 2 A and V Q = 32 V. This same set of probe- 
current measurements was also used as input to the program, but 
with the probe angle, <j» p , displaced ±5°. The computer output 
for these three sets of input data is presented in Figure 7-12, 
showing (reasonably so) that the principal error introduced by an 
uncertainty or error in the probe pitch-angle is in the thrust- 
loss factor, Ff Figure 7-13 shows the results of Figure 7-12 
plotted versus the error in probe angle, <*>p, and a linear 
approximation of these results gives the numerical value of the 
sensitivity as 3F . / 3 4> = 0.0017 deg -1 . 

In a similar manner, we evaluated the sensitivity of the 
thrust-loss factors to the radial coordinate, r, by changing the 
input values of the measurement coordinates. Figure 7-14 shows 
the output of the computer program (using the same input data 
described above), indicating that the principal effect of an 
error in the radial coordinate is manifest in the calculated 
value of the ratio J ++ /J + . (Fortunately, the thrust-correction 
factor a is relatively insensitive to this ratio.) The results 
of Figure 7-14 are plotted in Figure 7-15, and linearization 
gives the sensitivity as 3(J ++ /J+)/3r » 0.0096 cm" 1 . 

2. Uncertainty Estimates 

We estimated the uncertainties in the coordinates of 
Table 7-1 using the experimental setups indicated in Figures 7-16 
to 7-18. In the optical measurements depicted in Figure 7-16 we 
used a laser to calibrate the indicated probe pitch angle with 
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Figure 7-12. Computer printout showing the sensitivity of thrust-loss 
factors to an uncertainty in the angular coordinate, <b». 
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Figure 7-12. Continued 
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Figure 7-14. Continued 
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Figure 7-14. Concluded 
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the coordinate H. Then we calculated the angle <j> 0 using the 
geometrical relationship given by 

4»p = ~2~ tan-1 L(H - 92.7)/1 54.9 J . (7-16) 

The results of this calibration established the uncertainty in 
the angular coordinate as U^p = 0.45°. 

The uncertainty in the location of the probe with respect to 
the vacuum-chamber axis was determined using the laser and bore- 
siqht arrangement illustrated in Figure 7-17. The collimating 
tube established the location of the tank axis, using the same 
positioning pins that are used to align the thruster-mounting 
flange during installation of a thruster, with the chamber axis 
determined in this manner, the laser beam was adjusted to allow 
it to pass through the small— diameter (0.3 cm) apertures in the 
collimator. We observed that the probe-positioning rod (probe 
translational axis) was illuminated by the laser light, establish- 
ing the uncertainty (conservatively) as U. » 0.6 cm - the radius 

P 

of the positioning rod. 

The uncertainty in the radial coordinate r was determined . 
using the precision dial indicator shown in the setup of Figure 
7-18. Using this arrangement and a stepping-motor calibration 
factor (that was established at the time of the initial installa- 
tion 7-1 of the probe in the vacuum chamber) of 514 steps = 3.625 
cm, we indexed the probe up and down a fixed number of steps and 
recorded the change in dial-indicator reading. These measure- 
ments showed that there is negligible hysteresis in the posi- 
tioning mechanism, h least-squares fit of our data established a 
calibration factor of 141.8 steps/cm, in excellent agreement with 
the original calibration (141.8 steps/cm x 3.625 cm = 513.99 
steps). The largest deviation of the dial-indicator reading from 
the calibrated movement established the uncertainty in the radial 
coordinate as U r ■ 0. 'o54 cm. 
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3. Overall Uncertainty 

A summary of the uncertainties in the variables of Table 7-1 
and the sensitivities of the thrust-loss factors a and Ft is 
presented in Table 7-2. These results establish the overall 
uncertainties as U a ■ 0.1% and UF t * 1.4%. Substituting 
these values of the uncertainties into Equation 7-15 gives the 
overall uncertainty in the thrust as 1.4% (assuming zero error in 
the calculated thrust). 


Table 7-2. 


Summary of Sensitivity/Uncertainty Measurements 

14617-1 


3a 

3F t 

3Xj 

3x. ( 



1.69-* cm -1 

3.6 -4 cm -1 

1.69“ 3 cm -1 

3.6~ 4 cm -1 

0 

0 



4.0“ 4 deg -1 

9.89 -3 deg -1 

-4.9“® deg -1 

1.73 -3 deg" 1 

-4.9 -5 deg -1 

1.73 -3 deg -1 

4.0“ 4 deg -1 

9.89~ 3 deg -1 

~Q 

—0 


N/A 
0.6 cm 


0.079 cm 


SUM 

4.33 x 10“ 7 

1.90 x 10 -4 

y/SUM 

6.58 x 10“ 4 

1.38 x 10~ 2 










SECTION 8 


CONCLUSIONS 

The ring-cusp discharge chamber equipped with a three-grid 
ion-extraction assembly has demonstrated highly efficient and 
stable operating characteristics over a wide range of beam 
current and net-to-total accelerating-voltage ratio. Specific 
operating characteristics that were demonstrated include: 

* Beam-current throttling over the range 1 A < < 5 A. 

» Baseline specific discharge power of = 125 eV/ion. 

» Cathode-flow-dependent plasma conductivity, which 

enables stable control of discharge voltage over the 
range 28 V < < 44 V by regulating cathode flow 

rate. 

® The ability to operate the three-grid ion-extraction 
assembly over a beam-current range of 1 A < J5 < 5 A 
and net-to-total accelerating-voltage ratios in the range 
0.2 < R < 0.8. 

Plasma-property measurements obtained with the ring-cusp 
magnetic-confinement configuration indicate that the boundary of 
the ion-production region in this discharge chamber is defined by 
the shape of the confining magnetic-field lines, as it is in 
other more-conventional geometries. Steep gradients in plasma 
properties are observed in the vicinity of the magnetic-field 
lines that pass near the cathode. Essentially no ionization 
occurs in the volume exterior to these "critical" field lines, 
and their shape is thought to be an important feature in 
achieving the high level of performance demonstrated with this 
magnetic-confinement arrangement. 

Two additional discharge-chamber designs that employed 
single-cusp and multipole boundary-magnetic-field arrangements 
demonstrated the advantages of this confinement scheme in 
producing a uniform plasma-density distribution, and their 


i-Ri:ci:niNr. page plank not filmed 


177 


m, 



INJENIIQNAUI BLANK 






i 

f 

r ' 

L 




performance was comparable to that of the highly optimized 
J-series design. 

Perveance measurements obtained with an advanced ion- 
extraction assembly demonstrated an improvement in dimensional 
stability (over that of the state-of-the-art design) can be 
achieved by using flexible mounts for supporting both the screen 
and accelerator electrodes. 

A high-emission-current cathode design was derived from 
scaling the J-series cathode for long-life operation at emission 
currents up to Jg = 40 A, and three 20,000 A-hr endurance tests 
of the scaled cathode were successfully completed. The results 
of insert-temperature and orifice-plate-erosion measurements 
have demonstrated that long-life cathode operation should be 
achievable using a flow-restricting aperture diameter of 
d Q > 1.27 mm. 

Measurements of the minimum breakdown voltage for 
individual sections of the state-of-the-art J-series propellant- 
flow electrical isolator show that the conservative design value 
of the section breakdown rating is Vb = 340 V. The breakdown 
voltage of multiple sections was found to vary linearly, so that 
the total voltage-standoff rating can be expressed as the product 
of the number of isolating sections and the section breakdown 
rating of 340 V. 

Estimates of the uncertainties involved in evaluating the 
thrust-loss factors « and Ft indicate that the overall error in 
the results should be less than 1.4%. 
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APPENDIX A 


PLASMA PROPERTIES 

We conducted Langmuir-probe measurements in the ring-cusp 
discharge chamber, while operating at cathode-emission currents 
that placed us on either side of the transition from net- ion to 
net-electron current collection by the screen electrode. An 
initial survey of the discharge chamber allowed us to determine 
the locations of gradients in the plasma, and we used these 
results to define the locations where probe data would be nec- 
essary for defining the extent of the ion-production region. 

Figure A-1 shows a grid indicating the locations that were 
surveyed, with the "dots" indicating the locations where the 
probe data were recorded and analyzed. The close spacing of the 
"dots" indicates the regions where gradients in the plasma prop- 
erties were observed, and (as the results presented in Section 2 
showed) it also indicates where the probe characteristics were 
best represented by a Maxwellian distribution function. 

In the central region where the plasma was relatively 
uniform (and, in particular, near the cathode) we generally found 
that the probe characteristics were best approximated by a 
current-voltage relationship that represents the collection of 
both monoenergetic and Maxwellian electrons (two-group 
assumption). Figure A-2 shows probe data plotted in semi- 
logarithmic form for three locations on the chamber axis, ranging 
from near the cathode (location 1-1 ) to near the screen electrode 
(location A-1). When probe measurements are presented in this 
form, the presence of a monoenergetic electron group is manifest 
as curvature of the the low-voltage part of the characteristic. 

As the figure indicates, the characteristics change from a curved 
one (two-group) to a straight-line (Maxwellian), with increasing 
distance from the cathode. Numerical analysis of a characteris- 
tic similar to the one labelled 1-1 in Figure A-1 (corresponding 
to location L-3), shows that a good fit is obtained under the 
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Figure A-1 . Locations that can be accessed by the Langmuir 
probe in the ring-cusp thruster. The mesh 
points marked by the symbols are the locations 
where probe data were recorded and analyzed . 
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two-group assumption, but a single-group (Maxwellian) analysis 
produced a poor fit. Figures A- 3 and A-4 present a comparison of 
the input current-voltage characteristic with the results 
obtained by the curve fit, confirming the validity of the two- 
group model (for this location, which is near the source of 
emitted electrons). 

In most of the region characterized by gradients in the 
plasma, the numerical-analysis algorithm was unable to obtain a 
good fit under the assumption of a two-group plasma. To ensure 
completeness, we analyzed the probe data using both the two-group 
and single-group assumptions. The results of the two-group anal- 
ysis are presented in Table A-1 , and the results of the single- 
group analysis are presented in Table A-2 . The plasma-property 
measurements were performed while operating the thruster at = 

2 A, J D ■ 8 A, and V_ = 36 V, (e, ■ 144 eV/ion) , corresponding to 
net-ion current collection by the screen grid. Results 
corresponding to net-electron current collection by the screen 
grid (Jg = 9 A, = 162 eV/ion) are presented in Table A-3 
(two-group analysis) and Table A-4 (single-group analysis). 
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Figure A-3. 


A comparison of probe data (curve) with 
calculated results (symbols) that were 
obtained using a two-group model to analyze 
the probe characteristic/ 
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Figure A-4 . 


A comparison of probe data (curve) with 

«h^ a 5 ed F esults . (symbols) that were 

anali^ d fhf ln9 i si £9 le “9roup model to 
analyze the probe characteristic. 




Table A-1 . Plasma Properties of the Ring-Cusp Thruster Operated 
with Net-Ion Current Collection by the Screen Grid 

(Two-Group Analysis) 


TEMPERATURE 


CURVE 

MAXWELLIAN 

PLASMA 

PRIMARY 

MAXWELLIAN 

PRIMARY MAXWELLIAN 

NP/NM 

ION 

ION 


FIT 

SLOPE 

POTENTIAL 

ELECTRON 

SATURATION 

ELECTRON ELECTRON 


DENSITY 

FLUX 





ENERGY 

CURRENT 

DENSITY DENSITY 





(EV) 

(EV) 

(V) 

(EV) 

(MA) 

(PER CC) (PER CC) 


(PER CC) (MA 

PER CC) 

TRACE NO. 

6.89 

5.90 

37.17 

30.98 

76.26 

0 . 1422E+1 1 0. 2145E+12 

0.066 

0 . 2287E+12 

6.369 

AI/250B 

6.16 

6.16 

37.54 

33.13 

79.11 

0 . 1053E*11 0 . 2177E+12 

0.048 

0 . 2282E+12 

6.444 

A3/250B 

6.80 

5.80 

37.28 

29.56 

74.47 

0.9161E-10 0.2112E+12 

0.043 

0 . 2204E+12 

6.024 

A5/250B 

6.08 

5.08 

36.78 

34.07 

62.84 

0.3037E*10 0.1904E-12 

0.016 

0.1934£*12 

4.882 

A7/1008 

4.34 

4.34 

36.18 

29.72 

51.34 

0.2616F.*10 0.1684E-12 

0.016 

0 . 1 710E* 12 

3.988 

A8/100B 

3.24 

3.24 

35.86 

16.97 

29.56 

0 . 5557E*10 0 . 1122E+12 

0.050 

0 . 1 177E+12 

2.411 

A9/100B 

3.40 

3.40 

36.16 

25.33 

35.38 

0 . 2015E+10 0.1311E-12 

0.015 

0 . 133 IE* 12 

2.748 

A9/50B 

6.42 

6.42 

38.52 

37.38 

93.62 

0.7237E+10 0 . 2523E*12 

0.029 

0 . 2596E+ 12 

7.414 

C3/250B 

5.84 

5.83 

38.13 

31.80 

88.90 

0 . 7656E+10 0.2514E-12 

0.030 

0 . 2591E+ 12 

7.058 

C5/250B 

5.66 

5.65 

38.31 

32.53 

86.94 

0.5211E-10 0 . 2497E*12 

0.021 

0 . 2550E+12 

6.804 

C9/250B 

3.36 

3.37 

35.93 

22.31 

56.29 

0 . 9280E+10 0. 2093E+12 

0.044 

0 . 2186E+12 

4.552 

C10/100 

5.70 

5.71 

38.65 

29.51 

96.61 

0.1088E*11 0 . 2762E*12 

0.039 

0 . 2871E+12 

7.768 

E5/250B 

3.51 

3.51 

37.18 

20.41 

63.53 

0 . 9654E+ 10 0 . 2316E* 12 

0.042 

0 . 2413E+12 

5.123 

E8/2S0B 

6.74 

6.75 

39.81 

35.89 

113.66 

0 . 2036E+1 1 0 . 2988E+12 

0.068 

0 . 3191E+12 

9.518 

G1/250B 

5.99 

5.99 

39.42 

33.81 

104.86 

0.7886E+10 0 . 2927E+12 

0.027 

0. 3006E+12 

8.281 

G3/250B 

4.25 

4.25 

37.90 

25.41 

83.31 

0 . 1045E»11 0 . 2761fc '12 

0.038 

0 . 2865E+12 

6.684 

G8/250B 

8.38 

8.39 

41.13 

36.19 

126. 0 V 9 

0.2704E+11 0 . 2973E+12 

0.091 

0.3243E-12 

10.895 

I1/250B 

3.55 

3.54 

37.31 

22.22 

80.82 

0.1765F>11 0.2934E-12 

0.060 

0 . 3110E+12 

6.699 

I3/250B 

9.48 

9.51 

42.05 

35.12 

132.80 

0. 5074E+11 0 . 2941E*12 

0.173 

0.3448E+12 

12.774 

K1/250B 

5.14 

5.14 

38.47 

28.83 

85.83 

0 . 9658E+10 0.2586E+12 

0.037 

0 . 2682E-* 12 

6.881 

K2/250B 

9.53 

9.56 

41.28 

32.63 

108.39 

0.7106E*11 0.239SE+12 

0.297 

0 . 3106E+12 

12.134 

Ll/2503 

4.51 

4.51 

38.38 

23.67 

68.39 

0 . 1066E+11 0 . 2198E-*12 

0.049 

0 . 2305E«12 

5.571 

L2/250B 


T aoie 


with Net--i? n ^cS?Jen^colLction i brth? P s ThrUSter 0perated 

X clon °Y the Screen Grid 
( S ingle “Group Analysis) 


TEMPERATURE 
CURVE MAXWELLIAN 
p IT SLOPE 


(EV) 
3.22 
1.63 
8.01 
7.84 
7.09 
6.06 
5.26 
4.24 
4.46 
1.86 
8.11 
7 . 74 
6.92 
$ . 73 
6.60 
6.49 
5.37 
1.96 
8.36 
7. 57 
6.94 
6.01 
4.72 
2.84 
1.89 

9.21 

7.22 

5.70 
6.83 

5.71 


PLASMA 

POTENTIAL 


t?£!$EL WS»&fcJ*N 


(EV) 

(V) 

ENERGY 

(EV) 

CURRENT 

(MA) 

3.22 

36.90 

0.00 

17.07 

1.63 

36.93 

0.00 

2.03 

8.01 

38.70 

0.00 

89.62 

7.84 

38.67 

0.00 

89.20 

7.09 

38.15 

0.00 

82.76 

6.06 

37.88 

0.00 

66.34 

5 . 26 

37.36 

0.00 

54.34 

4.24 

36.82 

0.00 

33.99 

4.46 

37.72 

0.00 

37.43 

1.86 

37.10 

0.00 

9.91 


8.11 

7.74 

6.92 

6.73 

6.60 

6.49 

6.37 

1.96 

8.36 

7.57 

6.94 

6.01 

4.72 

2.84 

1.89 

9.21 

7.22 

5.70 
6.83 

5.71 


ma **ELLIAN 

DENSITY 
(PER CC) (PER CC) 


NP/NM 


ION 

DENSITY 


ION 

FLUX 


39.75 

39.54 

38.89 

39.15 

39.16 
38.97 

38.57 

37.17 

40.58 
39.94 
39.54 
39.08 
38,41 
37.35 
37.24 
41.47 

40.62 

39.62 
40.20 
39.42 


0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


103.02 
101 .22 
96.10 
97.52 
97.16 
91.99 
66.34 
14.18 
Ho. 49 
110.67 
106.40 
94.99 
71.63 
35.19 
15.93 
133.81 
113.09 
96.41 
109.59 
93.03 


C.0000E+00 0.6497E+11 
0. 0000E+00 0 . 1084E* 1 1 
0.0000E+00 0,2163E*12 
0.0000E+00 0.2176E*12 
0 . 0000E-*00 0.2122E*12 
0 . 0000E+00 0 . 1840E+12 
0.0000E+00 0. 1618E*12 
0.0000E+00 0.H27E+12 
0.0000E+00 0. 1210E+12 
0.0000E+00 0 . 4968E+11 
0.0000E*00 0 . 2471E+ 1 2 
0 . 0000E+00 0.2484£*12 
0 . 0000E+00 0.2495E+12 
0 . 0000E+00 0.2568E+12 
0 . 0000E+00 0.2583E+12 
0 . 0000E+00 0.2466E+12 
0 . 0000E+00 0 . 1926E+12 
0.0000E+00 0.6925E-»ll 
0.0000E400 0, 2752E+12 
0 . 0000E+00 0.2748E412 

0.0000E+00 0.2758E412 
0 . 0000E+00 0. 26465*12 
0.0000E400 0.2252E412 
0.0000E400 0 . 1427E+12 
0 . 0000E+00 0.7923E411 
0.0000E+00 0. 3012E+12 
0.0000E400 0.2876E412 
0.0000Ef00 0.281BE412 
0.000OE400 0.2863E412 
0.0000E400 0.2669E*12 


0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.000 


0 . 1084E* 1 1 
0 . 2163E* 12 
0 • 21 76E*12 
0 . 2122E+12 
0. 1840E+12 
0.161 8E* 1 2 
0 . 1 12?E*12 
0. 1210E+12 
0 . 4968E+1 1 
0 • 24715*12 
0.24845*12 
0 . 2495E + 12 
0 . 2568E+12 
0 . 2583E+12 
0.2466E412 
0.1926E412 
0.692SE411 
0. 2752E+12 
0.2748E412 
0. 2758E+12 
0. 2646E+12 
0.22S2E412 
0.14275*12 
0.7923E4H 
0.3012E412 
0.2876E412 
0.281SE412 
0.2863E412 
0.2659E412 


(MA PER CC) 

trace NO 

1 1.296 

A10/600 

0.154 

All/108 

6.801 

A1/260B 

6.769 

A3/2S06 

6.280 

A6/250B 

5 .038 

A7/1008 

4.124 

A8/1000 

2.579 

A9/1008 

2.841 

A9/50B 

0.752 

C11/S0B 

7.819 

C1/2S0B 

' .682 

C3/250B 

7.293 

CS/250B 

7.401 

C7/250B 

7.373 

C8/2S0B 

6.982 

C9/250B 

4.969 

C10/100 

1.076 

E10/50B 

3.840 

El/2508 

8.399 

E3/250B 

8.076 

E5/2S08 

7.209 

E7/2S0B 

5.436 

E8/250B 

2.670 

E9/100B 

1.209 

C10/S0B 

10.155 

C1/250B 

8.683 

O3/260B 

7.468 

C5/250B 

8.317 

07/2608 

7.060 

08/2608 
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Table A- 2 . Continued 


TEMPERATURE 
CURVE MAXWELLIAN 

PLASMA 

PRIMARY 

MAXWELLIAN 

PRIMARY 

MAXWELLIAN 

NP/NM 

ION 

ION 

FIT 

SLOPE 

potential 

ELECTRON 

SATURATION 

ELECTRON 

ELECTRON 

DENSITY 

FLUX 

cev) 

(EO 

(V) 

ENERGY 

(EV) 

CURRENT 

(MA) 

density 
(PER CC) 

DENSITY 
(PER CC) 


(PER CC) (MA 

, PER CC) 

3.21 

3.21 

37.92 

0.00 

43,49 

0 . 0000E*00 

0 . 16586*12 

0.000 

0. 16586*12 

3.300 

9.70 

9.70 

40.54 

0.00 

146.93 

0.0000E+00 

0.32226*12 

0.000 

0 . 32226*12 

11.151 

7.48 

7.48 

40.49 

0.00 

116.67 

0.0000E+00 

0 . 29136*12 

0.000 

0.29136*12 

8.854 

5.61 

5.61 

39.47 

0.00 

96.27 

0 .00006*00 

0 . 2776E+12 

0.000 

0 . 2776E+12 

7.306 

4.02 

4.02 

38.63 

0.00 

74.43 

0.0000E+00 

0.2535E+12 

0.000 

0 . 2536E*12 

5.649 

3.19 

3.19 

38.28 

0.00 

55.78 

0.0000E+00 

0.2132E*12 

0.000 

0 .21326*12 

4.233 

2.70 

2.70 

37.94 

0.00 

48.31 

0.0000E+00 

0 . 2007E-* 12 

0.000 

0.20076*12 

3.667 

2.42 

2.42 

37.72 

0.00 

41.93 

0 .0000E+00 

0.18416*12 

0.000 

0 . 1841E*12 

3.182 

2.37 

2.37 

37.78 

0.00 

39.42 

0 . 0000E+00 

0 . 1748E+ 12 

0.000 

0.1748E*12 

2.992 

2.16 

2.16 

37.80 

0.00 

32.42 

0 .0000E+00 

0 . 1507E+12 

0.000 

0.15076*12 

2.461 

9.43 

9.43 

38.51 

0.00 

166.40 

0 . 0000E*00 

0 . 3700E+12 

0.000 

0.37006*12 

12.628 

9.43 

9.43 

38.51 

0.00 

166.40 

0 . 0000E+00 

0 . 3700E+12 

0.000 

0 . 3700E*12 

12.628 

9.43 

9.43 

38.51 

0.00 

166.40 

0 .0000E+00 

0 . 3700E*12 

0.000 

0.37006*12 

12.628 

6.53 

6.53 

39.75 

0.00 

94.97 

0.0000E-00 

0 . 2537E*12 

0.000 

0.25376*12 

7.207 

6.53 

6.53 

39.75 

0.00 

94.97 

0.0000E+00 

0 . 2537E*12 

0.000 

0 . 2537E*12 

7.207 

6.53 

6.53 

39.75 

0.00 

94.97 

0.0000E+00 

0 . 2537E*12 

0.000 

0 . 2537F.*12 

7.207 

3.99 

3.99 

38.65 

0.00 

57.73 

0 . 0000E+00 

0.1974E*12 

0.000 

0.19746*12 

4.381 

3.99 

3.99 

38.65 

0.00 

57.73 

0.00006*00 

0 . 1974E+12 

0.000 

0 . 1974E*12 

4.381 

3.99 

3.99 

38.65 

0.00 

57.73 

0.0000E+00 

0 . 1974E*12 

0.000 

0.19746*12 

4.381 

2.70 

2.70 

37.86 

0.00 

41.04 

0 . 00006*00 

0 . 1705E*12 

0.000 

0.17056*12 

3.114 

2.70 

2.70 

37.86 

0.00 

41.04 

0.00006*00 

0 . 1705E*12 

0.000 

0. 17056*12 

3.114 

2.70 

2.70 

37.86 

0.00 

41.04 

0.0000E+00 

0 . 1705E*12 

0.000 

0 . 1705E*12 

3.114 

1.95 

1.95 

37.33 

0.00 

29.25 

0.00006*00 

0.1432E*12 

0.000 

0.14326*12 

2.220 

1.95 

1.95 

37.33 

0.00 

29.25 

0 . 00006*00 

0 . 1432E+12 

0.000 

0 . 1432E*12 

2.220 

1.95 

1.95 

37.33 

0.00 

29.25 

0 . 00006*00 

0 . 1432E*12 

0.000 

0.14326*12 

2.220 

1.57 

1.57 

37.61 

0.00 

21.06 

0.0000E+00 

0. 1147E*12 

0.000 

0. 1147E*12 

1.598 

1.29 

1.29 

37.31 

0.00 

15.78 

0.00006*00 

0,94836*11 

0.000 

0.94836*11 

1.198 

1.12 

1.12 

37.16 

0.00 

12.06 

0 .0000E+00 

0.78026*11 

0.000 

0.76026*11 

0.916 

1.03 

1.03 

37.18 

0.00 

10.34 

0 . 00006*00 

0.69676*11 

0.000 

0.69676*11 

0.784 

9.14 

9.14 

36.90 

0.00 

154.23 

0. 00006*00 

0.3484E*12 

0.000 

0.34846*12 

11.705 
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trace no. 

G9/100B 

11/2508 

12/2608 

13/2606 

I4/250B 

16/1008 

I6/100B 

17/1008 

18/1008 

I9/50B 

K1/250B 

K1/250B 

K1/250B 

K2/250B 

K2/250B 

K2/250B 

K3/100B 

K3/100B 

K3/100B 

K4/100B 

K4/100B 

K4/100B 

K5/50B 

KE/50B 

K5/50B 

K6/50B 

K7/50B 

K8/50B 

K9/50B 

L1/250B 



Table A- 2 . Concluded 


TEMPERATURE 


CURVE 

MAXWELLIAN 

PLASMA 

PRIMARY 

MAXWELLIAN 

PRIMARY 

MAXWELLIAN 

NP/NM 

ION 

ION 


FIT 

SLOPE 

potential 

ELECTRON 

SATURATION 

ELECTRON 

ELECTRON 

density 

FLUX 





ENERGY 

CURRENT 

OENSITY 

OENSITY 





(EV) 

(EV) 

(V) 

(EV) 

(MA) 

(PER CC) 

(PER CC) 


(PER CC) (VA 

PER CC) 

TRACE NO 

6.72 

5.72 

39.28 

0.00 

77.22 

0 . 0000E+00 

0 . 2205E + 12 

0.000 

0 . 2205E* 12 

5.860 

L2/2S0B 

3.32 

3.32 

36.39 

0.00 

43.01 

0.0000E+00 

0 . 1813E+ 12 

0.000 

0 . 1613E+12 

3.264 

L3/100B 

2.15 

2.15 

37.65 

0.00 

28.60 

0.0000E+00 

0 . 1333E+12 

0.000 

0 . 1333E+12 

2.170 

L4/100B 

1.65 

1.65 

37.68 

0.00 

19.75 

0.0000E*00 

0 . 1050E+12 

0.000 

0 . 1050E+12 

1.499 

L5/60B 

1.16 

1.16 

37,02 

0.00 

13.38 

0 . 00005*00 

0 . 8480E+11 

0.000 

0 . 8460E+11 

1.016 

L6/50B 

1.10 

1.10 

37.13 

0.00 

10.88 

0.0000E+00 

0 . 7076E+ 1 1 

0.000 

0 . 7076E+1 1 

0.826 

L7/50B 


Table a- 3 . Plasma Properties of the Ring-Cusp Thruster Operated 
with Net-Electron Current Collection by the Screen Grid 

(Two-Group Analysis) 


TEMPERATURE 
CURVE MAXWELLIAN 
FIT SLOPE 

(EV) (EV) 

PLASMA 

POTENTIAL 

(V) 

PRIMARY 

ELECTRON 

ENERGY 

(EV) 

MAXWELLIAN 

SATURATION 

CURRENT 

(MA) 

PRIMARY MAXWELLIAN 

ELECTRON ELECTRON 
DENSITY DENSITY 

(PER CC) (PER CC) 

NP/NM 

ION 

density 

(PER CC) (MA 

ION 

FLUX 

PER CC) 

TRACE NO. 

8.66 

8.66 

38.56 

28.05 

82.43 

0.7278E-10 0 . 1914E+12 

0.03e 

0. 1966E+12 

6.616 

A1/260C 

8.57 

8.67 

38.41 

33.49 

83.79 

0 . 2851E+10 0. 1956E+12 

0.015 

0.1983E+12 

6.499 

A3/250C 

6.18 

6.18 

36.96 

28.44 

58.77 

0. 5011E+10 0 . 1615E+12 

0.031 

0. 1665E+12 

4.669 

A7/250C 

8.81 

8.81 

39.18 

28.65 

95.59 

0 . 3409E+10 0.2200E+12 

0.015 

0 . 2234E+12 

7.42* 

C1/250C 

7.94 

7.94 

36.90 

32.41 

90.12 

0 . 707SE+10 0 . 2184E+12 

0.032 

0 . 2254E+12 

7.170 

C3/250C 

7 .04 

7.04 

38.35 

30.46 

83.45 

0.9829E*10 0 . 2148E+12 

0.046 

0 . 2246E+12 

6 . / 72 

C5/250C 

6.86 

6.86 

38.37 

32.17 

85.63 

0.6973E-10 0. 2233E+12 

0.031 

0 . 2303E*12 

6.807 

C7/250C 

7.10 

7.10 

38.62 

31.95 

88.77 

0 . 4807E+10 0.2275E*12 

0.021 

0.2323E-12 

6.950 

C8/250C 

6.78 

6.78 

38.41 

34.33 

84.23 

0 . 447 1E+ 10 0 . 2209E+12 

0.020 

0 . 22546+12 

6.586 

C9/250C 

7.97 

7.97 

39.65 

35.54 

101.26 

0. 1105E+11 0.2450E*12 

0.045 

0 . 2561E *12 

8.215 

E1/2S0C 

7.76 

7.76 

39.41 

33.26 

98.92 

0 . 7760E+10 0 . 2425E+12 

0.032 

0.25035*12 

7.871 

E3/250C 

7.61 

7.61 

39.29 

36.15 

97.72 

0 . 2996E+10 0 . 2420E+ 12 

0.012 

0 . 2450E+12 

7.558 

E5/250C 

6.38 

6.38 

38.74 

31.23 

84.97 

0.7426E*10 0 . 2298E+12 

0.032 

0 . 2372E+12 

6.766 

E7/250C 

6.33 

5.32 

38.34 

30.18 

68.21 

0 . 3761E+10 0. 2019E+12 

0.019 

0 . 20576*12 

5.325 

E8/100C 

8.19 

8.21 

40.29 

35.37 

110.83 

0 . 1921E+11 0.2642E+12 

0.073 

0 . 2834E* 12 

9,337 

G1/250C 

7.69 

7 59 

39.71 

35.47 

105.84 

0 . 3944E+ 10 0 . 2625E+12 

0.015 

0 . 26646*12 

8.213 

G3/250C 

5.31 

5.31 

38.47 

28.03 

83.98 

0 . 1162E+11 0 . 2489E+12 

0.047 

0 . 2605E+12 

6.824 

G5/250C 

6.79 

6.79 

39.46 

32.24 

96.95 

0.8364E+10 0 . 2542E+12 

0.033 

0 . 2625E+12 

7.723 

G7/250C 

5.46 

5.46 

38.18 

30.44 

83.01 

0 . 9247E+10 0 . 2426E+12 

0.038 

0 . 25196*12 

6.662 

G8/250C 

2.94 

2.95 

36.96 

18.29 

38.30 

0 .6117E+10 0.1524E412 

0.040 

0 . 1585E+12 

3.082 

G9/100C 

9.72 

9.74 

41.67 

35.16 

124.66 

0 . 2949E-*11 0.2728E*12 

0.108 

0 . 3023E*12 

11,024 

I1/260C 

8.03 

8.02 

40.12 

30.29 

110.58 

0 .6809E+10 0 . 2666E*12 

0.026 

0.2734E*12 

8.716 

I2/250C 

6.44 

5.43 

38.60 

30.56 

86.10 

0 .8777E+10 0 .2523E+12 

0.035 

0.261lE*12 

6.879 

I3/250C 

3.80 

3.80 

38.26 

17.87 

55.07 

0.2176E+10 0 . 1929E+12 

0.011 

0 . 19S0E+12 

4.247 

I5/100C 

9.12 

S . 13 

40.83 

33.49 

107.43 

0 . 7686E+11 0 . 2428E+12 

0.317 

0.3197E*12 

12.306 

K1/500C 

7.27 

7.27 

39.68 

36.48 

93.08 

0 . 1170E+10 0.2358E+12 

0.005 

0 .2370E*12 

7.118 

K2/2S0C 

3.86 

3.86 

37.54 

25.81 

56.56 

0 . 4634E*10 0 . 1967E+12 

0.024 

0 . 20136*12 

4 . 448 

K3/250C 

9.72 

9.73 

40.10 

31.73 

91.67 

0.9820E+11 0.2007E+12 

0.489 

0 .2989E*12 

12.637 

L1/500C 

6.30 

6.31 

38.95 

31.93 

78.14 

0 . 5598E+10 0 . 2125E+12 

0.026 

0.2181E+12 

6.163 

L2/250C 

3.73 

3.73 

38.06 

16.36 

41.63 

0.4831E+10 0. 1472E+12 

0.033 

0 . 1520E+12 

3.313 

L3/100C 


191 


Table A-4 . Plasma Properties of the Ring-Cusp Thruster Operated 
with Net-Electron Current Collection by the Screen Grid 

(Single-Group Analysis) 


temperature 


CURVE 

FIT 

(EV) 

MAXWELLIAN 

SLOPE 

(EV) 

PLASMA 

POTENTIAL 

(V) 

PRIMARY 

electron 

ENERGY 

(EV) 

4.25 

4.25 

37.05 

0.00 

2.07 

2.07 

37.40 

0.00 

9.06 

9.06 

38.61 

0 .0b 

8.92 

8.92 

38.58 

0.00 

8.20 

8.20 

38.04 

0.00 

7.02 

7.02 

37.58 

0.00 

6.44 

6.44 

37.24 

0.00 

5.57 

5.57 

37.05 

0.00 

6.46 

6.46 

38.75 

0.00 

2.37 

2.37 

37.15 

0.00 

8.84 

8.84 

38.97 

0.00 

8.50 

8.50 

38.96 

0.00 

8.11 

8.11 

38.91 

0.00 

7.79 

7.79 

38.96 

0.00 

7.68 

7.68 

39.16 

0.00 

7.60 

7.60 

39.12 

0.00 

2.72 

2.72 

37.48 

0.00 

9.20 

9.20 

40.22 

0.00 

8.51 

8.51 

39.70 

0.00 

8.04 

8.04 

39.55 

0.00 

7.38 

7.38 

39.43 

0.00 

6.17 

6.17 

39.20 

0.00 

3.72 

3.72 

37.76 

0.00 

2.56 

2.C6 

37.44 

0.00 

9.51 

9.51 

40.28 

0.00 

8.11 

8.11 

40.09 

0.00 

6.80 

6.80 

39.68 

0.00 

7.78 

7.78 

40.09 

0.00 

7.04 

7.04 

39.78 

0.00 

4.20 

4.20 

38.43 

0.00 


MAXWELL IAN PRIMARY MAXWELLIAN 
SATURATION ELECTRON ELECTRON 

CURRENT DENSITY DENSITY 


(MA) 

(PER CC) (PER CC) 

15.80 

0 , 0000E-00 0.5232E-11 

1.55 

0 , 0000E-00 0.7342E-10 

88.18 

0.0000E-00 0 . 2001E-12 

86.35 

0.0000E+00 0 . 1975E-12 

80.19 

0 .0000E+00 0. 1913E-12 

63.29 

0.0000E-00 0.1632E-12 

53.31 

0 .0000E-00 0 . 1435E-12 

34.70 

0.0000E*00 0.100SE-12 

63.47 

0. 00006-00 0 . 1 705E-12 

10.21 

0.0000E-00 0 . 4530E-11 

97.95 

0.0000E-00 0 . 2251E-12 

95.96 

0.0000E-00 0 . 2249E-12 

91.99 

0. 00006-00 0 . 2206E-12 

92.03 

0 . 0000E-00 0 . 2253E-12 

93.08 

0.0000E-00 0.2294E-12 

88.78 

0 . 0000E-00 0 . 2200E-12 

14.57 

0.0000E-00 0.6034E-11 

111.35 

0.0000E-00 0, 2508E-12 

105.67 

0.0000E-00 0 . 2475E-12 

100.60 

n . 0000E-00 0.2424E-12 

91.81 

0 . 0000E-00 0.2308E-12 

71.91 

0 . 0000E-00 0 . 1978E-12 

34.48 

0.0000E-00 0 . 1221E-12 

16.20 

0 . 0000E-00 0.6918E-11 

126.66 

0.0000E-00 0 . 2806E-12 

109.71 

0 . 0000E-00 0 . 2630E-12 

94.50 

0.0000E-00 0. 2475E-12 

104.58 

0 . 0000E-00 0.256 IE- 12 

92.36 

0.0000E-00 0. 2378E-12 

43.48 

0.0000E-00 0 . 1449E-12 


NP/NM 

ION 

density 

ION 

FLUX 



(PER CC) (MA PER CC) 

TRACE NO. 

0.000 

0 . 5232E-11 

1.199 

A10/50C 

0.000 

0.7342E-10 

0.117 

A11/10C 

0.000 

0. 2001E+12 

6.692 

A1/250C 

0.000 

0 . 1975E-12 

6.554 

A3/250C 

0.000 

0 . 1913E-12 

6.086 

A5/250C 

0.000 

0 . 1632E-12 

4.803 

A7/250C 

0.000 

0.143SE-12 

4.046 

A8/250C 

0.000 

0 . 100SE- 12 

2.633 

A9/100C 

0.000 

0 . 1705E-12 

4.817 

C10/100 

0.000 

0.4530E-11 

0.775 

CU/50C 

0.000 

0 . 2251E-12 

7.434 

C1/250C 

0.000 

0.2249E-12 

7.282 

C3/2S0C 

0.000 

0 . 2206E-12 

6.981 

C5/250C 

0.000 

0 . 2253E-12 

6.984 

C7/250C 

0.000 

0. 2294E-12 

7.064 

C8/250C 

0.000 

0 . 2200E-12 

6.737 

C9/250C 

0.000 

0.6034E-11 

1.105 

E10/50C 

0.000 

0 . 2508E-12 

8.451 

E1/250C 

0.000 

0.247SE-12 

8.020 

E3/2S0C 

0.000 

0.2424E-12 

'.634 

E5/250C 

0.000 

0.2308E-12 

6.968 

E7/250C 

0.000 

0 . 1978E-12 

5.458 

E8/100C 

0.000 

0.1221E-12 

2.617 

E9/100C 

0.000 

0 .69186-11 

1.230 

G10/S0C 

0.000 

0. 2806E-12 

9.612 

G1/2S0C 

0.000 

0.2630E-12 

8.326 

G3/260C 

0.000 

0.2475E-12 

7.172 

G5/250C 

0.000 

0 . 2561E-12 

7.937 

G7/2S0C 

0.000 

0.2378E-12 

7.008 

G8/250C 

0.000 

0 . 1449E-12 

3.300 

G9/100C 
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Table A-4 


Continued 


TEMPERATURE 


CURVE 

FIT 

(6V) 

MAXWELLIAN 

SLOPE 

(EV) 

PLASMA 

POTENTIAL 

FRIMARY 

ELECTRON 

ENERGY 

(EV) 

MAXWELLIAN 

SATURATION 

CURRENT 

(MA) 

FRIMARY 
ELECTRON 
DENSITY 
(PER CC) 

MAXWELLIAN 
ELECTRON 
DENSITY 
(PER CC) 

NP/NM 

ION ION 

DENSITY FLUX 

(PER CC) (MA PER CC) 

2.34 

2.34 

37.67 

0.00 

25.17 

0 .0000E+00 

0 . 1123E+12 

0.000 

0 . 1123E+12 

1.910 

1.29 

1.29 

37.06 

0.00 

10.01 

0 . 0000E+00 

0 .6015E+11 

0.000 

0.60156+11 

0.760 

1.29 

1.29 

37.06 

0.00 

10.01 

0.0000E+00 

0.60156+11 

0.000 

0.6015E+11 

0.760 

9.69 

9.69 

39.08 

0.00 

142.61 

0.0000E+00 

0 .31296+12 

0.000 

0.31296+12 

10.823 

8.36 

8.38 

40.11 

0.00 

1x5.86 

0 . 0000E+00 

0 . 2734 E+ 12 

0.000 

0 . 2734E+12 

8.793 

6.79 

6.79 

39.82 

0.00 

94.64 

0 . 0000E+00 

0.2480E+12 

0.000 

0 . 2480E+12 

7.182 

5.11 

5.11 

39.01 

0.00 

73.92 

0 . 0000E+00 

0 . 2233E+12 

0.000 

0 . 22336+12 

5.610 

4.01 

4.01 

38.44 

0.00 

56.70 

0 . 0000E+00 

0.1934E+12 

0 . 000 

0. 19346*12 

4.303 

3.24 

3.24 

38.10 

0.00 

46.36 

0 .0000E+00 

0 . 1760E+12 

0.000 

0 . 1760E+12 

3.518 

2.96 

2.96 

37.90 

0.00 

40.78 

0.0000E*00 

0.1618E+12 

0 . 000 

0 . 16186+12 

3.095 

2.79 

2.79 

37.81 

0.00 

37.61 

0 . 0000E+00 

0 . 1538E+12 

0.000 

0 . 15386+12 

2.854 

2.62 

2.62 

37.83 

0.00 

33.28 

0 . 0000E + 00 

0.14056+12 

0.000 

0.14056+12 

2.526 

0.85 

0.85 

37.12 

0.00 

5.63 

0 . 00006+00 

0.4185E+11 

0.000 

0.41856+11 

0.428 

9.26 

9.26 

37.40 

0.00 

167.91 

0 . 0000E+00 

0 . 37706+12 

0.000 

0. 3770E+12 

12.743 

9.26 

9.26 

37.40 

0.00 

167.91 

0 . 0000E+00 

0 . 3770E+12 

0.000 

0.3770E+12 

12.743 

11.18 

11.18 

38.94 

0.00 

162.90 

0 . 0000E+00 

0.3328E+12 

0.000 

0.33286+12 

12.363 

11.18 

11.18 

38.94 

0.00 

162.90 

0 . 0000E+00 

0 . 332BE+12 

0.000 

0. 3328E+12 

12.363 

7.45 

7.45 

39.69 

0.00 

94.21 

0. 0000E+00 

0.2358E+12 

0.000 

0.2358E+12 

7.150 

4.81 

4.81 

38.52 

0.00 

60.92 

0 . 0000E+00 

0 . 1897E+12 

0.000 

0 . 1897E+12 

4.623 

3.20 

3.20 

37.84 

0.00 

40.02 

0.0000E+00 

0. 1528E+12 

0.000 

0.15286+12 

3.037 

2.39 

2.39 

37.53 

0.00 

28.23 

0 . 0000E+00 

0 . 1248E+12 

0.000 

0 . 1248E+12 

2.142 

1.91 

1.91 

37.40 

0.00 

20.35 

0.00006+00 

0. 1006E+12 

0.000 

0.10066+12 

1.544 

1.61 

1.61 

37.49 

0.00 

14.42 

0.0000E+00 

0,77646+11 

0.000 

0.7764E+11 

1.094 

1.17 

1.17 

36.87 

0.00 

10.38 

0.0000E+00 

0.6549E+11 

0.000 

0.65496411 

0.788 

0.84 

0.84 

37.16 

0.00 

4.71 

0.00006*00 

0.3504E+U 

0.000 

0.35046+11 

0.358 

0.84 

0.84 

37.15 

0.00 

4.71 

0.0000E+00 

0 . 35046+11 

0.000 

0.3S04E+U 

0.358 

9.08 

9.08 

36.62 

0.00 

173.74 

0.0000E*00 

0. 3939E+12 

0.000 

0.39396412 

13.186 

10.38 

10.38 

36.35 

0.00 

156.98 

0.0000E+00 

0.3329S+12 

0.000 

0.33296412 

11.914 

7.24 

7.24 

39.69 

0.00 

83.53 

0 .0000E+00 

0 . 2121E+12 

0.000 

0.21216412 

6.339 

4.15 

4.15 

38.31 

0.00 

45.10 

0.0000E+00 

0 . 1513E+12 

0.000 

0.15136412 

3.423 


TRACE NO. 
I10/50C 
I11/S0C 
I11/50C 
I1/2S0C 
I2/250C 
I3/250C 
I4/250C 
IS/100C 
I6/100C 
I7/100C 
I8/100C 
I9/50C 
K10/10C 
K1/250C 
K1/250C 
K1/500C 
K1/500C 
K2/2S0C 
K3/2S0C 
K4/100C 
K5/50C 
K6/50C 
K7/50C 
K8/50C 
L10/10C 
L10/10C 
L1/2S0C 
L1/600C 
L2/250C 
L3/100C 
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Table ^-4 


Concluded 


TEMPERATURE 


CURVE 

FIT 

MAXWELLIAN 

SLOPE 

PLASMA 

POTENTIAL 

PRIMARY 

electron 

ENERGY 

MAXWELLIAN 

SATURATION 

CURRENT 

PRIMARY 

ELECTRON 

DENSITY 

MAXWELLIAN 

ELECTRON 

DENSITY 

NP/NM 

ION 

DENSITY 

ION 

flux 

TRACE NO 

(EV) 

(EV) 

(V) 

(EV) 

(MA) 

(PER CC) 

(PER CC) 


(PER CC) (MA PER CC) 

2.60 

2.60 

37.61 

0.00 

23.06 

0 . 0000E+00 

0.1189E+12 

0.000 

0. 1189E+12 

2.130 

L4/100C 

2.60 

2.60 

37.61 

0.00 

28.06 

0.0000E+00 

0 . 1189E+12 

0.000 

0 . 1189E+12 

2.130 

L4/100C 

1.87 

1.67 

37.47 

0.00 

18.71 

0.0000E+00 

0 . 9347E*11 

0.000 

0 . 9347E+11 

1.420 

LS/50C 

1.87 

1.87 

37.47 

0.00 

18.71 

0 .0000E+00 

0 . 9347E+1 1 

0.000 

0.9347E-11 

1.420 

L6/50C 

1.41 

1.41 

37.19 

0.00 

12. S2 

0.0000E+00 

0 . 7205E+11 

0.000 

0 . 7205E+11 

0.950 

L6/50C 

1.41 

1.41 

37.19 

0.00 

12.52 

0 . 0000E+00 

0 . 7205E+11 

0.000 

0 . 7205E+11 

0 . 950 

L6/50C 

1.21 

1.21 

37.06 

0.00 

9.41 

0 .000CE+00 

0 . 5831E+11 

0.000 

0 . 5831E+1 1 

0.714 

L7/50C 

1.21 

1.21 

37.06 

0.00 

9.41 

0.0000E+00 

0.5831E+1! 

0.000 

0. S831E+11 

0.714 

L7/50C 

1.06 

1.06 

36.83 

0.00 

8.67 

0.0000E+00 

0 . 5752E+ 1 1 

0.000 

0 . 5752E+1 1 

0.658 

L8/50C 

1.06 

1.06 

36.83 

0.00 

8.67 

0 .0000E+00 

0 . 5752E+11 

0.000 

0 . 5752E+11 

0.658 

L8/50C 

0.91 

0.91 

36.76 

0.00 

7.90 

0 . 0000E + 00 

0 . 5664E+11 

0.000 

0.5664E-11 

0.600 

L9/50C 

0.91 

0.91 

36.76 

0.00 

7.90 

0 , 0000E+00 

0 . 5664E+1 1 

0.000 

0 . 5664E+11 

0.600 

L9/50C 
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APPENDIX B 


DERIVATION OF PERVEANCE CORRECTION FACTOR 

The soace-charqe-limited current J+ that can be extracted 
by one-dimensional planar electrodes is governed by Child's 
Law: 


J 


+ 





(3-1 ) 


where V>j» is the total voltage, l e is the effective inter- 
electrode spacing, e/m is the charge-to-mass ratio, A s is the 
effective area, and e Q is the permittivity. Equation (3-1) 
applies to a beam composed of singly charged ions; if all the 
beam ions were doubly charged, the current would fc^ greater by 
a factor of JT. Therefore, in making porveance comparisons, 
the increase in measured perveance caused by the presence of 
doubly charged ions should be accounted for by correcting the 
measured value to obtain an equivalent singly charged ion 
perveance. 

For the case where both singly and doubly charged ions are 
present, Poisson's equation is 


dx 2 


JL 

e 


(J + + 


m 

2e 


V 


- 1/2 


(B-2 ) 


which can be integrated to give the appropriate form of Child's 
Law describing space-charge-limited flow of a beam composed of 
singly and doubly charged ions: 


e 


i 


y . 
f ■ 




‘.p 


r 





Rearranqing Equation (B-3) qives the familiar expression for 
the pe.r-ve^mc.e. .of a planar diode 


j + + j ++ /yr 

~ W 2 


2 
e 

The measured perveance is given by 



(B-4) 


( B-5 ) 


_ _ _± ++ ,o_ fi X 

P measured “ ,,3/ 2 ' 

V T 

where the sum of the singly and doubly charged ion currents is 
the measured beam current Jfc>. Combining Equations (B-4) and 
(B-6) gives the following relationship between the effective (or 
corrected) perveance and the measured value: 

P „ , J . + J.. A/2 

yorreg^d „ ±_ -±± v . (B -7) 

measured + ++ 

The right side of this expression is the definition of the 
thrust correction factor, a, so that Equation (B-7) can be 
rewritten as 


P 


corrected 


a P 


measured 


(B-8 ) 


APPENDIX C 



DERIVATION OF EXPRESSIONS FOR THE THRUST-LOSS FACTORS a AND F fc 


Expressions for the thrust-loss factors a and F^ can be derived by 
applying the integral form of the steady-state continuity and momentum 
equations to the control volume of Figure C-l. The continuity equation 
for the flow of ions* across the control surface can be written as 

J b = J + + °++ 

= f (J + + J ++ ) • ndA , (C-l) 

where J. i s the beam current, 1 is the current-density vector, n is the 

g 

unit vector normal to the small surface element dA, and the integration 
is over the electrode boundary A^ of the control volume. The scalar 
product appearing in Eq. C-l can be written as 

j • n ■ j cos (y - 4>) , (C-2) 

where * and y are the inclination of the current-density vector and 
surface normal with respect to the thruster axis. The angle 4 is 
determined experimentally; the angle y is determined by the electrode 
geometry and is given by the expression 

y = sin ^ ( jp j , (C-3) 

c 


where r is the radial coordinate, and R Q is the radius of curvature of the 
electrodes. The area element dA is related to the radial coordinate r 
(radial symmetry is assumed) by the expression 


* We assume that the ion beam contains singly and doubly charged 
ions only. 



iitifiti 


’sntar.muOi 


b) GEOMETRY OF SURFACE-AREA E LEMENT d A 


Figure C. 


1 . Control volume and coordinates used in 
deriving expressions for beam current 
and- net thrust. 
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•At * 



(C-4) 


dA = 


2irrdr 
cos y ’ 


which is obtained by projecting the area dA onto a circular plane 
perpendicular to the thruster axis. The thrust is obtained from the 
momentum equation as 



where m is the mass of the propellant ion, e is the electronic charge, 
and v + is the velocity of singly charged ions. In formulating Eq. C-5, 
we have neglected pressure forces as well as thrust produced by the loss 
of neutral particles and the impingement of charge-exchange ions. The 
factor 2 ~ in Eq. C,-5 accounts for the higher velocity and lower mass-to- 
charge ratio of doubly charged ions. The scalar product in Eq. C-5 
can be written as 


v + • n = v + cos ( Y - <j>) . (C-6) 

The component of thrust parallel to the thruster axis F is obtained 
by using the component j * z in Eq. c-5 

+ ir^++) ’ * cos (y - 4>)dA, (C-7) 


IK 

+ 2 “ -3++) cos <P cos (y - d>)dA . (C-8) 



199 


i*T?r.rTr> 


— n 


Alternatively, the net thrust is obtained by forming the scalar product 
of the thrust vector given by Eq. C-5 with the unit vector z 



where F t is the cosine of the average off-axis angle i. 

Equating the two expressions for F^ et results in an expression for the 
thrust-loss factor F t 


A„ 

F. = — 2- 


j (j + + J++) cos (y " cos 


f (J + + y- j ++ ) cos ( Y - 

ft 


(C-10) 


*) dA 


Combining Eqs. C-l and C-9 gives the net thrust in terms of the 
singly and doubly charged components of beam current J + and J ++ 


mv, 

Net " e 


- -r F t (J + + f J ++ > • 


(C-ll 


The calculated thrust is given by the expression below which neglects 
the losses resulting from off-axis trajectories and doubly charged ions 

mv + 

^Calculated i~~ J b 
mv. 

■ ~r (J + + J ++ > • (c-i2 

Net thrust is obtained by correcting the calculated thrust using the 
expression 


* 


Net 


t Calculated * 


vu- 1 


- ar 


Combining Eqs. C-ll, C-l 2 , and C-13 gives the following expression for 
the thrust-loss factor a 


a = 


J + + T J ++ 


J + + J H 


(C-V 


Eqs. C-l, C-l 2, and C-13 can be combined to give 


m 


Net ' aJ b e V + F t • 


(C-l! 


Substituting v 


A 


2eV 


b ... 


m 


gives the following expression for net thrust 


in terms of beam current and beam voltage 


F Net ■ -Wr \\ ■ 


(C-ll 


In this form, aJ K represents the equivalent singly charged ion current 
2 D 

and F t V. is the effective beam voltage. 



APPENDIX D 


DERIVATION OF EQUATIONS USED TO ANALYZE ExB PROBE DATA 


The probe collector current i corresponding to species n is given 
by the expression 

i n ( r > ? p ) =/ f I n (r,i) cos' - j{u[:-( J 0 -,)]-u[i-( y*)]}T(<p p -q;d w dA p , (D-l) 

where I n is *he intensity ;f "^e ion flux, r and ; are the coordinates, $p 
is the probe angle, v is the acceptance half-angle of the collimator, u is 
the unit step function, T is the transmission of a cylindrical collimator, 
a) is the solid angle, and A p is the aperture area. The combination of the 
transmission and unit step functions in Eq. D-l effectively "collimates" the 
incoming flux (as illustrated in Figure D-l), allowing only a fraction of 
those ions having angles in the range $ + $ to reach the velocity filter. 

In practice the operator selects the ion species of interest by varying the 
ratio E/B to match the particle speed. The narrow acceptance angle (2ip) 
of the collimator permits the intensity in Eq. D-l to be replaced by the value 
corresponding to the probe angle Since the collimator restricts 
the angle difference to <p - <j>p <_ \j> the cosine can be replaced by unity, 
and Eq. D-l can be written as 


1 n (r ’*P ) * (D-2) 

The narrow field-of-view of the collimator and the small size of the 
entrance aperture permits us to assume that the ion flux is homogeneous 
over the aperture area and isotropic within the field-of-view. This 
enables the intensity to be expressed as 


A similar expression is frequently used to describe the response 
of charged-particle analyzers employed on exploratory spacecraft (Ref. D-l). 
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where J is the current density and n is the solid angle subtended by the 
collimator. Reference to Figure d- 2 shows that the solid angle n is 
given by 


n 



(D-4) 


where the aperture area Ap is given by 


V 


.2 ,2 

z ip 


(D-5) 


where we have used the approximation tan \f> = The differential solid 
angle du is given by 

dA_ 

du = — p , (D-6) 

z 

- 2iWR . (D-7) 

jr 


For small angles R = zc and the expression above can be written as 


du - Znede . 


(D— 8) 


Combining Eqs. D-2, D-3, D-4, D-5 and D-8 f and performing the 
integration over dAp, gives 


2J n (r, < b)A n 

, n(r> * ) = niJ 


I 




eTde , 


(D-9) 


which can be solved for the current density J n 
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Figure D.2. Definition of variables used in analysis of 

the response of a cylindrical collimator of ; 

length 2 t. 


i 


» 


I 






(0-10) 





The transmission of the cylindrical collimator is given by 

T= 1 “lj'T 1 [ 1 - {^r) ] / + sin’ 1 (1JS.) | , (D-n) 

where, from Figure D-3, 

to . 1 S1n (0-12) 

2 r p 

Equations 0-11 and D-12 give the transmission of a cylindrical collimator 
inclined at an angle ♦with respect to a parallel beam. In our case, we 
are interested in the response of the collimator to a superposition of beams 
(the isotropic and homogeneous assumption) which are inclined at an angle e 
with respect to the probe axis (the angle e varies from 0 to v>). Therefore, 
the app-opriac* angle to use in the transmission function is e, and for small 


angles 

the equation above reduces to 

1+S £e 


where 

2 " r ' 

P 

— » i 

(D-13) 


r p ♦ 

» < . j • « 

(D-14) 


The transmission represents the fractional overlap of the circular areas 
of the entrance and exit apertures (as shown in Figure D-3) and, there- 
fore, represents the response of the collimator to a parallel beam. 

Combining Eqs. D-ll, D-13, and 0-14 gives 

T ” 1 -f|i L l - (!) j ' ♦ sin" 1 (tjj . (0-15) 


This expression is taken from Ref. D-2. An error in the equation 
appearing in the reference has been corrected. 
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A. COLLIMATOR ARRANGED PARALLEL TO BEAM (T * 1 ) 




B. COLLIMATOR INCLINED AT AN ANGLE 
0p WITH RESPECT TO BEAM (0<T< 1 ) 


Figure D.3. 


?hf 1 I? iti0n - 0f • variables used in determining 
the transmission of a cylindrical collimator, 


:riry JTTfcfr. 
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Evaluation of the integral appearing in Eq. D-10 is simplifed by 
making the substitution 




which gives 


J ET(c)dE * /J x - |[x 2 yiI7 + X sIiT’x] dx . 


The integration can be performed using standard techniques to evaluate 
the following integrals: 


/ xdx-i, 


/ dx - is- ■ 

J n 


x sin xdx = g- . 


Substituting these results into Eq. d-17 


gives 


/ o eTU,dc - * 2 [?-7(k+f)] 


(D-2 


Combining Eqs. D-10 and D-21 gives the relationship between the current 
density at the probe aperture and the collector current 



The relationship between the current density at the accel electrode 
and the collector current can be derived using the geometric variables 
defined in Figure D-4. The small area A Q of the accelerator electrode 
is viewed" from different angles $ , resulting in a current-density 
dispersion profile similar to the one shown. This profile represents 
the current density at the measurement plane due to ion flow from the 
area A Q . The total current from A Q is obtained by integrating this 
current density over the measurement plane 


I 


n 



n dA , 


(D-23) 


where the scalar product is given by 


3 • n = J 


cos 


The area element dA is given by the expression 


dA = - 


dAj_ 


cos 4 > q » 


(D-24) 


( D— 25) 


which is obtained by projecting the area element dA onto a plane normal 

to the direction defined by the angle <j> 0 corresponding to the peak of 

the current-density distribution. In this plane, the current density is 

assumed symmetrical about the radii defined by the angles d> . and <t> 

min T max 

of Figure D-4. Under this assumption, the area element dA^ is given by 


dAj_ = itRdR 


(D-26) 


* 

Note that the current density 0 is not necessarily the total current 
density at the measurement plane, since the collimator restricts the probe 
viewing area to A 0 . For this reason, the total current from the acceler- 
ator system cannot be found by integrating the current density at the 
measurement plane. 
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VOLUME SURFACE MEASUREMENT 
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Figure D , * . Definition of variables used in analysis 
of probe-collector current. 
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where the coordinate R is given by 


R = L sec tan ( t ? - t Q ) , (D-27) 

and 

dR - L sec <p Q sec 2 U p - . (D-28) 

Note that the area element dA^ is a half-ring element of length rrR and 
width dR, which is consistent with the symmetry assumption. 

The distance L from the measurement plane to the accelerator grid is 
given by 


L = L o + M 1 " cos » (D-29) 

where L Q is the distance from the probe to the accelerator electrode on 
the thruster axis, and R c is the radius of curvature of the electrodes 
(as shown, R Q is positive for electrodes dished outward). The angle 
y is given by 


Y * Sln '’(^-) • (0-30) 

Combining Eqs. b-22 through b-28 gives the desired expression for the 
current through the measurement plane 

sec 3 <f> 0 C * max 2 

‘n * r p J cos *P (♦p-» 0 )|tan (* p -+ 0 )|d* p . (0-31) 

Vin 

This is also the current through the area A Q of the accelerator grid, 
which can be written as 


Ip ■ J n COS (y - *) * 0 


(D-32) 
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where j is the current density at the accelerator electrode. Combining 
Eqs. 0-31 and D-32 gives the desired expression relating the current 
density at the electrodes to the probe current 


V r) 


4irL^sec^<|> 0 
A q A C0S Iy-4>) 


j maX i n (r^ p ) r -os$ p 3ec 2 (^ p -^ 0 )|tan(4> p -<j> 0 )|d4> p • (D-33) 

Vin 


An expression for the angle <f> can be derived by applying the momentum 
equation to the control volumes of Figure D-4. Using the measurement 
plane as the control surface, the net thrust is given by 


F net^ 



(D-34) 


Using the accelerator electrode as the control surface, the net thrust is 
given by 


F ( r ) = ^ [j + (r) + y- j ++ (r)] cos (y - <t>) cos * A Q (D-35) 

lie w ^ 


Equating the two expressions for net thrust and solving for the angle 4» 
gives 


4 >(r) = cos 


-1 


/[ 0 + (r,<t> p ) + J ++ (r,<fr p )] cos 2 4 > p dA 

[j + (r) + jr* J ++ (»')3 cos (y - +)A 0 


. (D-36) 


Substituting Eqs. 0-22, D-25 through 0-28, and D-33 into Eq. D-36 gives 
the desired expression for the angle <(> 

~ *The thrust equation is derived in Appc?.4ixC ; , (Eq. C-8). 
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